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I.  INTRODUCTION 


This  report  is  concerned  with  the  analysis  of  unsteady  two- 
dimensional  transonic  flow  in  channels.  Although  the  motivation 
behind  the  work  was  the  study  of  the  effects  of  oscillating  walls  on 
the  flow,  it  was  found  that  in  the  problems  studied  it  was  possible  to 
include  the  effects  of  oscillations  in  the  back  pressure,  interpreted 
as  pressure  variations  impressed  upon  the  flow  at  some  point  down¬ 
stream  of  the  shock  wave.  Hence  either  of,  or  any  mixture  of,  the 
two  means  of  causing  unsteadiness  in  the  flow  may  be  considered;  of 
particular  importance  is  the  ability  to  describe  the  kinds  of  wall  motion 
which  correspond  to  a  given  oscillation  in  back  pressure.  Because  both 
modes  of  introducing  unsteadiness  are  retained,  the  results  have  appli¬ 
cations  to  the  flow  in  inlets  and  to  the  flow  between  the  compressor 
blades  in  turbojet  engines  and  to  the  diffusers  of  ramjet  engines,  as 
long  as  the  flow  is  in  the  transonic  regime. 

Physically,  wall  oscillations  may  be  associated  with  either  blade 
flutter  or  the  displacement  of  the  flow  due  to  a  turbulent  boundary  layer 
or  a  separation  bubble;  the  oscillations  in  back  pressure  are  associated 
with  pressure  fluctuations  caused  by  combustion  instabilities  or  inter¬ 
actions  with  the  wakes  of  blades  or  flow  asymmetries,  etc.  ,  downstream 
of  the  channel  in  question.  When  the  relative  velocity  at  the  tip  of  a 
rotor  blade  exceeds  the  sonic  value,  the  possibility  of  supersonic  un¬ 
stalled  flutter  has  been  demonstrated  [  1,Z]  .  In  addition,  interaction 
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between  the  core  flow  and  the  boundary  layer  has  been  observed  in 
experiments  involving  transonic  flows  through  channels  with  separation 
[3,4];  apparently,  self-induced  oscillations  occurred  with  variations  in 
the  position  of  the  shock  wave  and  the  location  and  size  of  the  separation 
bubble.  Although  these  two  phenomena  are  examples  of  cases  where 
the  effective  boundary  of  the  inviscid  flow  field  moves  with  time,  it  is 
clear  that  changes  in  back  pressure  could  certainly  affect  the  inter¬ 
actions.  Finally,  in  a  ramjet  engine  (and  conceivably  in  a  turbojet 
engine)  combustion  instabilities  could  propagate  upstream  and  thus  be 
the  source  of  the  oscillations  in  pressure  [5] .  Indeed,  in  a  ramjet 
engine,  it  appears  possible  that  self- sustained  oscillations  could  be 
generated  as  a  result  of  the  interaction  between  the  diffuser  and  the 
combustion  chamber;  when  an  internal  shock  wave  occurs  in  the  diffu¬ 
ser,  such  oscillations  could  lead  to  this  shock  being  disgorged  and  thus 
to  diffuser  "buzz"  or  even  more  severe  problems. 

The  experimental  work  done  at  McDonnell  Douglas  (e.g.  ,  Refs. 

[4]  and  [5])  has  led  to  analyses  of  viscous -inviscid  interactions  in 
transonic  internal  flows,  both  unsteady  [6]  and  steady  [7],  but  with  un¬ 
separated  flow;  results  have  not  predicted  self-induced  oscillations  and 
so  the  more  difficult  separated- flow  problem  is  now  being  considered 
by  Liou  [8]  ,  among  others. 

The  work  described  herein  is  concerned  with  one  aspect  of  the 
interaction  problems  mentioned  above;  as  such,  it  is  closely  related 
to  the  experimental  work  done  on  diffuser  flows  carried  out  by  Sajben 


? 


et  al  [4,5]  .  Specifically,  it  has  to  do  with  the  prediction  of  the  inviscid 
transonic  flow  field  with  a  shock  wave  in  a  two-dimensional  channel, 
using  asymptotic  methods  of  analysis.  Previous  to  this  work,  in 
studies  of  unsteady  transonic  flow  in  channels  [9,  10,  11,  12]  involving 
asymptotic  methods,  the  unsteadiness  was  caused  by  pressure  oscilla¬ 
tions  downstream  of  the  shock  wave  such  that  time  dependence  occurred 
only  in  second-  and  higher-order  terms;  the  flow  was  steady  to  first 
order.  The  solutions  found  are  valid  for  inviscid  flows  in  symmetric 
and  asymmetric  channels  with  small-  and  large-amplitude  oscillations 
in  shock -wave  position;  in  all  cases,  the  walls  are  stationary.  It  can 
be  shown,  also,  that  a  transonic  flow  between  the  blades  of  a  two-dimen¬ 
sional  cascade  is  described  asymptotically  as  a  channel  flow.  Because 
the  previous  channel-flow  solutions  were  written  in  terms  of  arbitrary 
wall  shapes,  they  may  also  be  applied  to  a  two-dimensional  cascade  flow 
with  arbitrary  airfoil  shapes  under  conditions  where  the  blades  are 
stationary  and  oscillations  in  pressure  originate  downstream  of  the 
cascade.  In  the  present  work,  the  same  range  of  parameters  is  con¬ 
sidered  as  in  Refs.  [9]  to  [l2] ,  but  for  the  case  where  the  unsteadiness 
could  be  caused  by  oscillations  in  wall  shape  or  back  pressure  or  by 
any  combination  of  the  two. 

The  fact  that  the  problems  considered  all  involve  inviscid  flow 
fields  should  not  be  construed  as  being  restrictive  insofar  as  use  of  the 
results  are  concerned.  Because  the  solutions  are  written  in  terms  of 


arbitrary  wall  shapes,  the  effects  of  boundary  layers  may  be  accounted 
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for  by  adding  their  displacement  thickness  to  the  actual  wall  shape;  in 
the  range  of  parameters  considered,  the  boundary  layer  is  quasi-steady. 
In  practice,  the  calculation  is  accomplished  by  setting  up  a  computation¬ 
al  procedure  in  which  the  boundary-layer  displacement  thickness  and 
channel-flow  conditions  are  calculated  simultaneously.  Such  a  pro¬ 
cedure  has  been  used  by  Liou  [7]  who,  in  calculating  the  pressure  dis¬ 
tribution,  also  accounted  for  the  interaction  between  the  shock  wave 
and  a  turbulent  boundary  layer  using  asymptotic  techniques  developed 
by  Messiter  [13]  ;  comparison  between  experimentally  derived  and  com¬ 
puted  pressure  distributions  showed  excellent  agreement. 

In  the  next  section,  the  general  problem  is  formulated  and  the 
specific  problems  are  defined  in  terms  of  the  relative  orders  of  the 
parameters.  The  following  sections  contain  the  solutions  obtained 
during  the  one-year  period  of  the  contract.  In  the  final  section,  the 
results  are  summarized. 


II.  GENERAL  PROBLEM  FORMULATION 


The  problem  under  consideration  is  that  of  an  open  duct  placed  in 
a  uniform  transonic  flow.  Inside  the  duct  the  cross-section  area  varies 
with  distance  along  the  duct.  On  the  outside,  the  duct  walls  are,  for 
convenience,  assumed  parallel  to  the  incoming  uniform  flow.  The  lack 
of  any  cowl  effects  on  the  outer  flow  is  not  important  to  the  results 
obtained;  such  effects  are  omitted  only  to  make  calculations  as  simple 
as  possible.  A  sketch  of  the  duct  and  the  notation  employed  is  shown  in 
Fig.  1.  The  coordinates  x  and  y  are  made  dimensionless  with  respect 
to  L,  where  the  total  width  of  the  duct  is  taken  to  be  ZL  .  (Overbars 
denote  dimensional  quantities.  )  The  corresponding  velocity  components, 
u  and  v,  are  dimensionless  with  respect  to  the  critical  sound  speed  in 
the  incoming  flow,  a*.  Another  characteristic  dimension  of  the  channel 
is  its  length,  denoted  by  c.  The  channel  flow  can  be  associated  with 
that  between  the  blades  in  a  two-dimensional  cascade  with  zero  stagger 
angle;  the  length  of  the  channel  is  then  the  chord  of  the  blade  in  this 
equivalent  cascade  flow.  For  the  present  calculations,  the  flow  at  the 
exit  cross-section  might  be  considered  to  exhaust  into  a  plenum,  or  the 
duct  might  be  regarded  as  joined  smoothly  to  fixed  plane  walls  extending 
further  downstream. 

In  general,  the  nondimen sional  parameters  associated  with  the 
flow  unsteadiness  are  the  characteristic  time  (or  frequency)  and  ampli¬ 
tude  of  the  impressed  oscillations  in  pressure  or  in  boundary  position, 
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made  dimensionless  with  the  fundamental  parameters  of  the  basic  steady- 
channel  flow.  Thus,  if  T  ^  =  1/f  is  the  characteristic  time  associated 
with  walls  or  back  pressure  oscillating  at  a  frequency  f,  then  we  define 


T 


(1) 


Here,  L/ a  is  a  characteristic  flow  or  residence  time  in  the  channel, 
because  the  reference  velocity  for  this  transonic  flow  is  a*;  however, 

L/  a  *  is  also  the  order  of  the  time  taken  by  an  acoustic  signal  to  traverse 
the  channel.  The  parameter  r  may  be  related  to  the  familiar  reduced 
frequency  used  in  flutter  analysis,  by  noting  that  if  s  =  2L  =  cs  is  the 
blade  spacing  (s  is  dimensionless  with  respect  to  the  chord,  as  shown), 
then. 


T 


(2) 


where 


(2-irf)  ( c/2) 


—  * 
a 


(3) 


is  the  reduced  frequency  associated  with  the  blade  motion,  based  on  the 
half-chord. 

The  upper  and  lower  wall  shapes  are  written,  respectively,  as 


(yw)  =  1  -  £  f  (x)  +  a  G  (x,t) 

W  u  u 


(y  )  =-l+e  f  (x)  -  a  G .  (x ,  t ) 

w  |  If 

G  =  {(x  -  x  )  (3 (t)  +  G  (t)} 

u,l  c  p  u,  I 


(4a) 

(4b) 

(4c) 


where  f^(x)  id  f^{x'  ^e  the  steady-state  wall  shape,  with  f  =  f  =  0 
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at  the  leading  edges,  and  t  is  the  time  made  nondimensional  with  T 
The  motions  of  the  upper  and  lower  walls,  which  may  or  may  not  be  in 
phase  with  each  other,  are  denoted  by  G  and  G  ,  respectively;  x  is 
the  center  of  rotation,  with  a  |3(t)  being  the  instantaneous  angle  made 
by  the  outer  duct  wall  measured  relative  to  the  x  axis,  and  G^(t)  repre¬ 
sents  a  plunging  motion  independent  of  x.  The  thickness  of  the  "blades" 

2  2 

or  duct  walls  is  then  0{e  )  where  e  «  1;  another  way  of  interpreting 
€  is  to  note  that  the  nondimensional  radius  of  curvature  of  the  nozzle 
formed  by  the  duct  inner  wall  at  the  throat  is  0(e-^).  The  parameter  or, 
with  a  «  1 ,  orders  the  amplitude  of  the  motion  of  the  walls,  this 
depending  upon  the  problem  considered.  As  will  be  seen,  it  is  the 
relative  orders  of  or,  t,  and  €  which  determine  the  order  of  the  ampli¬ 
tude  of  the  shock-wave  motion  and  the  orders  of  the  velocity  and  pres¬ 
sure  perturbations  due  to  the  unsteadiness  in  the  flow.  Finally,  it  may 
be  noted  that  for  a  symmetric  channel  with  symmetric  wall  motions , 
fj(x)  =  fu(x)  and  G^x.t)  = 

Since  the  flow  is  transonic,  the  undisturbed-flow  Mach  number  M 

00 

differs  from  one  by  a  small  amount.  The  difference  M  -  1  is  taken  here 

00 

to  be  0(e ),  which  is  the  case  of  greatest  interest  because  the  steady-state 
variation  in  channel  width  is  O(e^)  and  so  the  changes  in  Mach  number  within 

the  channel  are  of  the  same  order  as  M  -1.  Since  a  is  taken  to  be  an  inte- 

00 

gral  power  of  6,  the  expansions  for  the  velocity  components  have  the  form 


,  2  3 

u=l+€Uj+e  u^  +  e  u^  +  •  •  • 

2  3 

V  =  €  Vj  +  €  Vj  +  •  •  • 


(5a) 

(5b) 
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where  Uj,  u^»  •  •  •  and  v^,  •  •  •  are  functions  of  x,  y,  and  t.  These 
representations  are  correct  throughout  most  of  the  channel,  except 
perhaps  for  thin  regions  close  to  the  leading  and  trailing  edges,  the 
throat,  and  the  shock  wave.  Corrections  for  such  regions  are  intro¬ 
duced  later  when  needed. 

If  the  governing  equations  are  written  in  dimensionless  form,  one 
might  choose  to  make  the  time  dimensionless  with  respect  to  the  char¬ 
acteristic  residence  time;  i.e.  ,  T  =  T/(L/  a  ).  However,  for  impressed 
oscillations  on  the  flow  T  =  0(T  ,  )  and  so,  as  noted  earlier,  we  intro- 

Cn 

duce  a  time  coordinate  t  made  nondimensional  with  T  ;  then  T  is 

cn 

related  to  t  by  a  coordinate  stretching  as  follows: 


=  T  -  rt 


(6a, b) 


where  t  is  defined  in  Eq.  (1). 

The  relations  for  the  shape  and  velocity  of  the  shock  wave  may  be 
written  as 

xs(T,y)  =  xgQ(t)  +  €Xj(t)  +  ...  (7a) 


dx  .  T  dx  n 

_ _ s_  _  _1_  sO 

Us  "  dT  T  [_  dt 


+  e 


dx  , 
s  1 

— —  +  •  .  . 

dt 


(7b) 


dt  ~  dt 
3/2 

That  is,  it  can  be  shown  [9]  that  9x  /9y  =  0(e  )  and  that  v  ,  the  com 

S  3 

ponent  of  the  shock  wave  velocity  in  the  y  direction,  is  of  high  enough 
order  that  it  may  be  neglected.  The  velocity  of  the  shock  wave  ug  is 
positive  in  the  positive  x  direction.  Because  the  fluid  velocity  is  near 
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*'T?— 


******  mm 


the  speed  of  sound  and  the  velocity  relative  to  the  shock  wave  must  be 

supersonic,  the  velocity  of  the  oscillating  shock  wave  must  be  small 

compared  to  the  sonic  velocity;  i.  e.  ,  |u  |  «  1.  Thus,  if  t  =  0(1), 

s 

for  example,  dxg<j/dt  =  0  so  that  the  amplitude  of  the  shock  motion 
cannot  be  0(1)  but  is  at  most  O(e),  from  Eqs.  (7). 

It  is  desired  to  find  solutions  for  the  velocity  components  valid  to 

2 

0(€  ).  This  allows  calculation  of  the  thermodynamic  variables  to  order 
2 

€  also.  Now,  as  will  be  seen,  it  is  necessary  to  calculate  the  term  of 
3 

order  e  in  v  and  to  apply  the  boundary  conditions  at  the  wall  in  order 

to  derive  an  equation  for  an  unknown  term  in  u^.  Hence,  governing 

equations  for  third-order  terms  must  be  derived.  Because  the  flow  is 

inviscid  and  starts  at  uniform  conditions,  it  is  isentropic  up  to  the 

shock  wave.  Because  the  flow  is  transonic,  shock  waves  are  weak,  and 

3 

the  jump  in  entropy  across  the  wave  is  0(g  );  however,  the  gradient 

in  entropy,  transverse  to  the  flow,  is  the  important  term  in  ordering 

the  vorticity,  and  this  gradient  in  the  y  direction  is  at  least  one  order 
4 

higher,  i.e.  ,  0(e  ),  downstream  of  the  shock  wave.  Therefore,  up  to 

3 

and  including  terms  0(e  ),  the  flow  is  ir rotational  upstream  and  down¬ 
stream  of  the  shock  wave,  so  that  a  velocity  potential  exists.  Thus, 
one  can  write  a  perturbation  potential  in  the  form 

2  3 

<}>  =  €  <(>j  +  e  +  e  +  •  •  •  (8) 

with  u  =  1  +  4>x,  v  =  <^,  Uj  =  4>lx  ,  Vj  =  etc.  The  governing  equa¬ 

tion  for  <$>  is,  then,  the  full  potential  equation  [14]  , 


(9) 


r 


2  2  2  2 

(a  -  u  )  <t>  +(a  -  v  )  <1>  -  2uv4  -  *_ 

XX  Tyy  Txy  tTT 


-  2  u  <j>xT  -  2  v  4>yT  =  0 


where  a  is  the  speed  of  sound  made  dimensionless  with  respect  to  a  *; 
a  is  given  by  the  Bernoulli  equation  for  unsteady  flow  [14]  : 


fn  +  4  +  *T  =  F,T| 


(10) 


2  2  2  2  2 

where  q  =  u  +  v  ,  and  a  /(V  -  1)  +  q  /2  =  H  is  the  stagnation  enthalpy, 


-*2 


dimensionless  with  respect  to  a  *  ;  F(T)  is  an  integration  function  to  be 
determined. 

A 

Expansions  for  the  density  p,  pressure  P,  and  temperature  T, 
each  made  dimensionless  with  respect  to  its  critical  value  (e.  g.  , 


p  =  p/  p  )  can  be  written  in  the  same  form  as  that  for  u,  Eq.  (5a);  thus. 


for  example. 


P  =  l+  €P1+e  p2  + 


(ID 


with  similar  expansions  for  P  and  T  .  Also,  it  is  easily  shown  that 


2  a 

a  =  T.  The  property  relationship,  in  dimensionless  form  is 


y  (IT)AS  Y  3 

p  =  py  e  =  pY  (1  +  0(e  )) 


(12a, b) 


where  Y  is  the  ratio  of  specific  heats,  As  is  the  difference  in  entropy 
* 


s  -  s  ,  made  dimensionless  with  respect  to  the  specific  gas  constant, 

R,  and  where  in  Eq.  (12b)  use  has  been  made  of  the  fact  that  at  most 
As  =  Ote3).  Thus,  from  Eqs,  (11)  and  (12b),  one  can  derive  the  relations 
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(13a, b) 
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u  are  dimensionless  with  respect  to  their  critical  values  in  the  incoming 
flow,  and  since  the  change  in  a  across  the  shock  wave  due  to  unsteady 
motion  is  at  most  of  order  e,  then  at  the  sonic  surface  and  u^  are 
both  zero  for  all  time,  and  the  result  therefore  becomes 


Pj  =  -  Uj 

Using  Eqs.  (19),  (13a),  and  the  first-order  potential  function,  one  can 


integrate  Eq.  (17b): 


H1  =  "  k  4* it  + 


where  is  a  function  of  integration.  If  Eqs.  (20)  and  (17a)  are  sub¬ 
stituted  into  the  definition  of  H  and  the  result  along  with  Eqs.  (8)  and 
(16)  is  substituted  into  Eq.  (10),  one  finds  that 


F(T>  -  2irrr> +  '2  «i + 


(21a) 


a2  =  1  +  (1  -q2)  +  (Y-  1)(H-H  )  (21b) 

2  oo 


and  g.  is  a  function  of  t  alone.  This  is  consistent  with  the  result  to  be 
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shown  later  that  4>j  is  independent  of  y.  The  potential  equation  (9)  can 

2 

now  be  rewritten  with  u  =  1  +  ,v  =  4>  ,  and  with  a  given  by  Eq.  (21b): 

x  y 


V  =  ml)  +  2  11  +  V+xT  *  *TT  + 


_  .  2+n  4. 

+  O  (e  ,  €  ) 


The  error  estimate  in  Eq.  (22)  uses  the  anticipated  results  that 

4>  =  O(e),  4>  =  0(e2). 
x  y 

Ahead  of  the  shock  wave,  following  a  streamline  upstream  to  the 
undisturbed  steady  flow,  one  finds  that  gj  =  0.  It  should  be  noted  for 
future  reference  that  there  is  no  change  expected  in  these  results  if  an 
inner  region  enclosing  the  leading  edges  of  the  duct  is  traversed  since 
no  new  oscillations  are  impressed  upon  the  flow  in  this  region. 

In  order  to  calculate  g^  downstream  of  the  shock  wave,  the  change 
in  H  across  a  moving  wave  must  first  be  calculated.  Since  this  change 
and  other  jump  conditions  are  required  for  shock  waves  moving  at 
velocities  of  varying  order,  it  is  necessary  to  consider  the  general 
jump  conditions  across  a  moving  wave.  These  conditions  are  derived 
most  conveniently  in  terms  of  a  coordinate  system  normal  and  parallel 
to  such  a  wave,  at  an  arbitrary  point  on  it,  as  shown  in  Fig.  2.  There, 
q  is  the  velocity,  n  and  t  are  unit  vectors  in  the  directions  normal  and 
perpendicular  to  the  wave,  respectively,  and  cn  is  the  normal  compo¬ 
nent  of  the  velocity  c  of  the  shock  wave  at  the  point  in  question;  the 
tangential  component  cfc  is  taken  to  be  zero.  If  the  superscript  +  is 
used  to  denote  conditions  relative  to  the  shock  wave,  and  double  brack¬ 
ets  [[  ]]  indicate  the  jump  in  a  quantity  across  the  wave,  then 

(qV  +  <qt+)2 

[[H+]]  =  [fh  +  —  -  —  3J  =  0  (23a) 


qt+  =  qt  (23b,  c) 
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2  ,  _  2 


Since  [[cn]]  =  0  and  H  =  h  +  (q^  +  q^.  )/2,  Eq.  (23)  can  be  written  as 


[fH]]  =  cn  [[q^]] 


and  it  remains  to  find  cr  and  q^. 


If  the  instantaneous  shape  of  the  shock  wave  is  denoted  by  x  (y,  T), 

s 

and  S  =  x  -  x  ,  then  S  =  0  on  the  shock  wave  and  also 
s 


S^,  +  c  •  VS  =  0 


(25a) 


-  _VS_ 

11  ~  l*S|  '  V  1  +  x  * 


(25b) 


where  i  and  j  are  unit  vectors  in  the  x  and  y  directions,  respectively, 
and  the  subscripts  y  and  T  denote  partial  derivatives.  Using  Eq.  (25b) 
for  VS  in  Eq.  (25a)  and  for  n  in  calculating  q^,  one  can  show  that 


c  = 
n 


(26a) 


1  +  x 


U  -  V  X, 


(26b) 


Also,  since  nXqsnXq^t, 


qt  = 


V  +  U  X 


1  +  X  2 

sy 


The  shock-wave  relations  for  one-dimensional  flow  are  applicable 
in  the  direction  normal  to  the  wave,  so  that  one  can  write  immediately, 


(qn+,d(qn+,u  =  (a+  ')Z 
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where  subscripts  u  and  d  refer  to  conditions  immediately  upstream  and 

+* 

downstream  of  the  wave,  respectively,  and  a  is  the  dimensionless 
critical  sonic  velocity  in  the  coordinate  system  fixed  to  the  wave.  Hence, 


(7  +  1) 
2(V-D 


1  .  ,  sT , 

- r-  x  (u  -  vx  )  -  — t—  ) 

2.  sT  sy  2 

1  +  x  7 

sy 


(29) 


3/2  2 

In  all  of  the  problems  considered  here,  xg^  =  0(e  ),  v  =  0(e  ),  u  =  0(1), 

and  x  _  =  O(e)  at  most.  If  the  shock-wave  velocity  in  the  x  direction  is 
si 

2 

defined  by  ug  =  xgT  .then  up  to  and  including  terms  0(e  )  one  can  write: 

c  ~  x  _  ~  u  (30a) 

n  sT  s 


u 


UX  +  V 

sy 


(u  -  u  )  ~  (a  ) 

s  u 


2(7-1) 


1  +*\2  U  .  3 

(a  )  ^  H  -  u  u  + 

3  & 


(30b, c) 
(30d) 

(30e) 


It  may  be  noted  that  xg  is  desired  only  up  to  and  including  terms  O(e), 
which  are  independent  of  y. 

If  the  inviscid  flow  field  is  considered  to  be  that  found  in  the  limit 

x 

as  the  viscosity  tends  to  zero,  such  that  the  continuity  of  <j>  =  f  d>  dx 

-00 

+  constant  through  the  shock  wave  is  retained,  then  4>  and  therefore  4>rp 
are  continuous  and  in  fact  unchanged  across  the  shock  wave.  Hence, 
F(T)  in  Eq,  (10)  (and  (21)  )  may  be  calculated  for  the  flow  downstream 
of  the  shock  wave  by  calculating  the  jump  in  H  across  the  shock  from 
Eq.  (24)  and  thus  the  jump  in  F(T)  using  Eq.  (10)  at  the  wave. 


The  boundary  conditions  to  be  used  at  each  order  of  approximation 


may  be  found  by  noting  that  if  W  =  y  -  yw(x,t),  then  at  the  surface  of  the 
duct  W  =  0  and  +  q  •  VW  =  0  there  also.  From  this  equation,  one 


finds  that 


v 

w 


u 

w 


dy 

1  w 

3x 


(31) 


where  the  subscript  w  denotes  conditions  at  the  wall,  or  surface  of  the 
duct. 

Finally,  specific  problems  are  defined  by  setting  the  relative 

orders  of  the  parameters  t ,  a,  and  e.  First,  we  choose  the  orders  of 

2  3 

magnitude  of  a  to  be  a  =  0(e  )  and  a  =  0(e  ).  That  is,  the  amplitude  of 
the  wall  motion  is  at  most  of  the  order  of  y^  -  1,  which  is  equivalent  to 
the  thickness  of  the  blades  in  a  cascade.  In  general,  the  order  of  mag¬ 
nitude  of  t  should  cover  a  range  from  t  » 1  to  t  =  0(1).  That  is,  for 
the  case  where  the  boundary  layer  apparently  interacts  with  the  moving 
shock  wave  the  flow  is  quasi-steady  so  that  t  »  1;  on  the  other  hand, 
when  the  unsteadiness  arises  because  of  blade  oscillation,  the  range  for 

r  includes  t  =  0(1).  For  the  problems  to  be  considered  here,  we  choose 
-2  -1 

t  =  0(e  )  and  t  =  0(€  ).  To  interpret  these  choices,  we  note  that  a 

weak  pressure  disturbance  travels  at  a  speed  |u-a|  ,  and  it  follows  from 

Eqs.  (5a),  (21),  and  (8)  that  u- a  =  O(c).  Since  the  distance  is  0(1), 

physically  the  case  v  -  0(e  *)  is  that  for  which  signals  from  the  source 

of  the  disturbance  reach  the  shock  wave  in  a  time  comparable  with  the 

time  characterizing  the  period  of  the  disturbance,  T  ,  ,  and  the  case 

ch 
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-2 

t  =  0(€  )  is  that  for  which  signals  reach  the  shock  wave  "instantaneously" 

compared  to  T  ,  ,  as  shown  in  Ref.  [  1 5 ]  .  The  case  t  =  0(1)  is  deferred 
ch 

to  future  work. 

The  specific  problems  considered  here  are  then: 

3 

(1)  Oscillations  of  Channel  Walls  with  a  =  0(e  )  and  of  Back  Pressure 
with  Perturbations  O(e^);  t  =  0(c-^). 

2 

In  this  case,  because  the  walls  are  stationary  to  0(e  )  and  the 

pressure  is  steady  to  O(e),  unsteadiness  in  the  flow  velocity  arises  only 

2  2, 

in  terms  0(e  ).  Hence  u  =  0(e  ),  and  from  Eq.  (7b)  it  is  seen,  there- 

s 

fore,  that  x  =  x  _  (t ) ;  i.  e.  ,  large -amplitude  motion  of  the  shock  wave 
so  s  0 

is  possible.  This  is  an  extension  of  the  problem  considered  in  Ref.  [12]. 

3 

(2)  Oscillations  of  Channel  Walls  with  a  =  0(e  )  and  of  Back  Pressure 
with  Perturbations  O(c^);  r  =  0(c"^) 

In  this  case,  the  flow  is  again  unsteady  in  second  order  but  the 

oscillations  occur  at  a  higher  frequency  than  in  problem  (1),  with  a 

2 

shock-wave  motion  of  smaller  amplitude;  i.  e.  ,  u  =  0(e  )  again,  so 

s 

from  Eq.  (7b)  xg^  =  constant.  Small -amplitude  shock  motion  results. 

2 

(3)  Oscillations  of  Channel  Walls  with  a  =  0(e  )  and  of  Back  Pressure 
with  Perturbations  0(e);  t  =  0(e“*). 

Here,  the  flow  is  unsteady  in  first  order,  and  so  u  =  O(c).  Hence, 

s 

as  in  problem  (1),  xg^  t  constant  so  large-amplitude  motion  of  the  shock 


wave  is  possible. 


Each  of  the  problems  listed  here  is  applicable  to  a  channel  with 


an  unsteady  boundary  layer  and/or  oscillating  back  pressure,  and  to 

blade  flutter,  although  for  the  last  application  problem  (1)  has  marginal 

3 

use.  In  problems  to  be  considered  in  the  future,  t  =  0(1)  with  a  -  (e  ) 

2 

in  one  case  and  a  =  0(e  )  in  the  other.  These  cases  will  be  applicable 
to  blade  flutter  with  marginal  application  to  unsteady  boundary  layers 
or  oscillating  back  pressures  because  of  the  small- amplitude  motion  of 
the  shock  wave. 
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in.  SOLUTIONS 


In  this  section  the  three  problems  described  in  section  II  are 

considered  in  detail.  In  general,  the  goal  has  been  to  obtain  solutions 
2 

valid  to  order  €  in  the  velocity  components  and  thus  in  the  pressure, 
temperature,  and  density.  As  the  problem  complexity  increases,  this 
is  not  possible  without  the  use  oi  numerical  methods  of  solution,  and  in 
those  cases  where  considerable  time  would  have  been  spent  in  obtaining 
these  higher-order  corrections,  the  problem  formulations  for  these 
terms  have  been  completed,  but  no  solutions  are  given. 

Problem  (1);  a  =  0(e  3),  t  =  0(c 

There  are  several  regions  in  which  the  main  channel-flow  solutions 
do  not  describe  the  flow  properly;  i.e.  ,  they  are  not  uniformly  valid 
there.  In  that  event,  special  "inner1*  solutions  which  hold  in  the  regions 
in  question  and  which  match  with  the  "outer"  channel-flow  solutions 
must  be  derived.  The  inner  and  outer  solutions  may  then  be  joined  to 
form  composite  solutions  valid  to  the  desired  order  of  magnitude  in 
larger  regions.  For  the  problems  considered  here,  inner  solutions 
are  needed  at  the  entrance  to  the  duct  (and  perhaps  at  the  exit),  at  the 
throat,  and  in  the  neighborhood  of  the  shock  wave.  In  addition,  a 
separate  calculation  must  be  made  for  the  velocity  and  position  of  the 
shock  wave.  Because  the  analyses  are  similar  in  all  of  the  problems 
to  be  considered,  each  of  these  special  calculations  is  contained  in  a 
titled  subsection  for  easy  reference. 


In  this  section,  the  equations  used  to  relate  a  to  e  and  r  to  f  are 


(32a) 


t  =  (k  e2)'1 


(32b) 


where  k  is  an  arbitrary  constant  of  order  one. 


Channel- Flow  Solutions 


In  order  to  find  the  governing  equations  for  each  order  of  approxi¬ 
mation  from  Eq.  (9),  it  is  necessary  first  to  calculate  a,  the  speed  of 

sound.  As  shown  previously,  F(T)  =  —  ( Y  +  1)/ ( Y  -  1 )  upstream  of  the 

2 

shock  wave  and  so  a  may  be  found  from  Eq.  (10)  immediately.  Down¬ 
stream  of  the  shock  wave  the  jump  in  H  is  found  using  Eqs.  (24),  (30a) 

2 

and  (30b),  with  u  =  ke  dx  /dt  +  •  •  •  from  Eq.  (7b);  thus,  one  can 
s  sU 

show  that  and  g^,  defined  in  Eq.  (21),  are 


a  dx 

H  -  H  =  e  k  [u.  .  -  u.  ] 

d  u  dt  Id  lu 


(33a) 


g!  =0 


(33b) 


Since  only  terms  0(e  )  are  needed  in  Eq.  (9),  only  terms  up  to  and 

2  2  1 
including  0(€  )  need  be  retained  in  a  ,  and  so  F(T)  =  —  (Y  +  1 )  /  (Y  -  1 ) 

suffices  throughout  the  channel  in  this  case. 

If  Eq.  (10)  is  substituted  into  Eq.  (9),  and  Eq.  (8)  substituted  into 

the  resulting  equation,  one  finds  the  following  governing  equations  for 

2  3 

the  4^  by  gathering  terms  of  order  e  ,  e  ,  and  e  ,  respectively,  and 
rearranging  slightly: 
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(34a) 


<t>,  =0 

lyy 

*3yy=  [(''tlH>lX*2x+  2k4>lt+ 


(34b) 

(34c) 


where  subscripts  x,  y,  and  t  denote  partial  derivatives.  If  Eqs.  (5), 

with  of  course  u  =  1  +  4>  and  v  =  <j>  ,  Eqs.  (4),  and  Eqs.  (32)  are  sub- 

x  y 

stituted  into  Eq.  (31),  giving  the  boundary  conditions,  one  finds  that  at 
the  upper  and  lower  walls,  respectively, 


2y  (*.l,t)  =  “  fu  * 

<|»3y(x,l,t)  =  GUx-u1f^ 

2y(x,-l,t)=  f;  , 

<)>  (x,  -l,t)  =  -  (G,  -  u 

3y  lx 

(35a, b) 

(35c, d) 


where  f 1  =  df  /dx  ,  etc. 
u  u 

From  Eq.  (34a),  since  there  is  no  term  O(e)  in  the  boundary  con¬ 
ditions  , 

=  4^  (x,t)  (36) 

Because  <j>^  is  independent  of  y,  Eq.  (34b)  may  be  integrated  easily. 
Upon  satisfying  the  boundary  conditions  on  4>2y,  given  by  Eqns.  (35a, c), 
one  finds  an  equation  involving  only  4^x>  f^,  and  f^  which,  although  non¬ 
linear,  again  is  easily  integrated  to  give: 


*lx  ±  V  (V+  1)  (C1  "  f) 


(37a) 


f  +  f. 

f  =  — L 

1  2 


(37b) 


where,  in  general,  Cj  =  Cj(t).  Now,  is  set  by  flow  conditions  at  the 
throat,  the  duct  exit,  or  the  duct  entrance.  For  example,  if  the  flow  is 
choked,  4>^  =  0  at  the  throat.  On  the  other  hand,  if  the  flow  throughout 
the  duct  is  subsonic,  conditions  at  the  exit  set  c^,  and  if  the  entrance 
flow  is  supersonic,  conditions  at  the  entrance  are  used.  In  the  case  of 
primary  interest  here,  the  entrance  flow  is  subsonic  and  choked  at  the 
throat  with  a  shock  wave  downstream  of  the  throat.  It  will  be  shown 
later  that  because  f  is  independent  of  time  and  only  second-order  pres¬ 
sure  oscillations  are  considered,  Cj  is  independent  of  time  for  this 
problem.  Hence,  if  we  set  x  =  0  at  the  point  where  f'  =  0, 


Cj  =  f(0) 


(38) 


Moreover,  <f>lt  =  0,  and  using  Eq.  (37a)  for  <j>lx,  the  solution  for  may 
be  written  simply  as 


2 

=  -  f*  +  g2  y  +  h2  (x,t) 


g2  = 


(f;  ‘  fu> 


(39a) 

(39b) 


where  h^  is  a  function  of  integration. 

With  cf>j  and  <\>^  known,  Eq.  (34c)  may  be  integrated  once  to  give 
<J>3^.  If  the  boundary  conditions  are  applied  at  y  =  ±  1  (Eqs.  (35b, d) ) 
and  the  resulting  equations  are  subtracted,  a  governing  equation  for 
h^(x,t)  is  obtained.  This  equation  may  be  integrated  once  to  give  h^, 
which  is  all  that  is  needed  for  the  velocity  components.  The  result  is 
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til 


(40a) 


c(t)+G 

h2*=  6  ['"-,2V-3)*1*]  +  — 


lx 


G  = 


G  +  G 

u  I 


(40b) 


In  Eq.  (40a),  c^(t)  is  a  function  of  integration  which  is  set  by  given  flow 

conditions  at  some  station  in  the  duct;  in  general  it  has  different  values 

upstream  and  downstream  of  the  shock  wave. 

With  these  results,  the  solutions  for  the  velocity  components  may 

2 

be  written  up  to  and  including  terms  0(e  ): 

,  2 

u=l  +  €Uj+e  u2  +  « « • 


=  1±£  (C~,  +  £2["  f”i‘  +  g2y  +  h2*)  +  " 


v=€  v2  + 


=  e 


[-£*  y  +  g7]  + 


(41a) 

:4  b) 


If  the  expansions  for  u  and  v  are  substituted  into  Eq.  (10)  with 

1  2 
F(T)  =  —  ("Y  +  1)/(V  —  1),  an  expansion  for  T  =  a  results ;  this ,  with 

Eqs.  (13a, b),  (14),  and  (19),  allows  the  following  relations  to  be 

derived  for  the  density,  temperature,  and  pressure: 

p=l-eu1-e2(u2+  u“)  +  ...  (42a) 

2 

A  2  Ul 

Tsl-efY-lJUj.  €  (V-l)(u2  +  -j-)  +  ...  (42b) 

2 

P=l-€\Uj-€  V  u2  +  •  •  •  (42c) 


Thus,  the  solutions  found  may  be  used  to  calculate  velocity,  demity, 
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temperature  and  pressure  distributions  across  the  channel  at  a  given 
location  in  the  flow  or  along  the  centerline  or  on  the  wall.  They  are 
valid  in  the  main  part  of  the  channel. 

It  is  seen  from  Eqs.  (42c),  (41a),  and  (40a),  that  imposing  a  given 
pressure  oscillation  of  order  at  the  exit  of  the  duct,  say,  where  x  =  x 
is  equivalent  to  prescribing  c^  +  G  there.  Thus  if  we  take  f  =  f*  =  f"  =  0 
at  x  =  Xj,  and  further  set 


Pz  (xrt} 


=  u2(Xj,t)  = 


ILL 


(s)  x 

3)  2  ,  .  C2d  + 

U1  (X1)  +  FTThT 


Gj(t) 


1(X1} 


(43) 


then  the  last  term  describes  the  oscillations  in  back  pressure  imposed 

(s) 

on  the  flow.  The  constant  c refers  to  a  steady-state  back  pressure, 

giving  a  steady- state  location  for  the  shock  wave,  and  G  ,(t)  refers  to  the 

d 

unsteady  fluctuations  in  back  pressure  superimposed  over  the  steady- 
state  value.  Then,  downstream  of  the  shock  wave,  c2(t)  is  set  by  the 
following  equation: 

c2(t)  +  G(Xj  ,t)  =  c2(®}  +  Gd(t)  (44) 


Upstream  of  the  shock  wave,  if  the  flow  is  choked, 

c2(t)  +  G(0,t)  =  0  (45) 

if  h2x  is  to  remain  finite  there  (see  Eqn.  (40a)  ).  This  is  discussed  in 
greater  detail  later,  and  it  is  also  shown  that  the  solution  for  u  can  be 
written  in  composite  form  such  that  the  singular  behavior  of  h2x  at  the 


throat  does  not  arise. 


It  is  important  to  note  that  the  oscillations  in  back  pressure  were 

Z 

chosen  here  to  give  flow  fluctuations  of  the  same  order,  namely  0(e  ), 

3 

as  those  induced  by  the  oscillations  of  the  wall  for  a  =  0(e  ).  This  was 
done  so  that  relevant  comparisons  of  the  two  effects  could  be  made. 

Leading-Edge  Solution 

Near  the  leading  edge  x  =  x^  the  "channel"  solutions  derived  above 
are  no  longer  correct  because  the  flow  is  no  longer  nearly  one-dimen¬ 
sional.  A  direct  matching  of  the  channel  solutions  with  the  corresponding 
undisturbed -flow  quantities  in  general  is  not  possible.  This  was  demon¬ 
strated  in  Ref.  [l6]  for  steady  flow  at  high  subsonic  speed  through  a 
cascade  with  staggered  blades  having  wedge-shaped  leading  edges,  with 
a  simple  correction  included  for  rounded  edges.  To  illustrate  the  need 
for  an  inner  leading-edge  solution  in  the  present  context,  the  simpler 
case  of  symmetric  cusped  edges  is  chosen,  so  that  f^fx^)  =  fj(xQ)  =  ® 
and  f"(xA)  =  f"(xj  =  f"(x.)  >  0;  wall  motions  are  such  that  G  (xn,t) 

=  G^(Xg,t)  =  G(Xg,t);  also  f^x^)  =  *j(xq)  =  0  by  definition.  The  leading- 
edge  solution  is  found  to  be  necessary  for  providing  correct  flow  details 
near  the  edges,  but  is  not  needed,  at  least  in  the  present  applications, 
for  determination  of  the  integration  constants  c^  and  c^  in  the  channel 
solution.  Solutions  for  wedge-shaped  edges  would  be  somewhat  more 
complicated,  and  would  lead  to  the  same  general  conclusion. 

For  symmetric  cusped  leading  edges,  expansion  of  the  solutions 
(41a, b)  for  u  and  v  as  x  —  x^  gives,  for  the  subsonic  case  where 
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<j>j' (x)  <  0  in  the  forward  part  of  the  duct. 


+  •  •  •  (46a) 

v  =  -  e2  f"  (xQ)  (x  -  xQ)  y  +  •  •  •  (46b) 

The  undisturbed-flow  velocity  is  represented  (Fig.  1)  as  a  series  expan¬ 
sion  in  e  by 

«oo=1+euloo  +  €2u2oo  +  ...  (47) 

where  u.  <0  for  the  subsonic  case.  The  velocity  perturbation 
X  oo 

u  -  1  =  4>x  is  also  0(e)  in  the  leading-edge  region,  and  the  term 

("Y  +  1)  <j>x  4>^  in  the  differential  equation  (22)  is  no  longer  small  in  com- 

1/2 

parison  with  <t>^  when  x  -  x^  =  0(€  ),  corresponding  to  a  thin  vertical 

strip  containing  the  edges.  For  flows  of  greatest  practical  interest, 

appreciable  velocity  changes  in  a  small  leading-edge  region  are  probably 

undesirable,  and  so  we  consider  only  cases  in  which  the  O(e)  term  is 

constant  there.  Comparison  of  equations  (46a)  and  (47)  then  allows  Cj 

to  be  related  to  u.  ,  which  is  considered  to  be  given: 
loo 


by 


jj(  j{c 

The  leading-edge  flow  depends  on  coordinates  x  and  y  ,  defined 


x  -  x. 
*  0 


<K€) 


1/2  * 


y  =  y 


(49) 
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where  the  factor  K  is  included  for  convenience,  with  K  defined  by 


1  -  M  2  =  Ke  +  0(e2)  =  -  (V  +  l)  u.  €  +  0(€2)  (50) 

00  loo 


The  perturbation  potential  has  the  form 


4>  =  (u  -  1)  x  +  e5^2  <j>  *  (x*.  y*,  t)  + 
00  4 


(51) 


A  complex  velocity  can  therefore  be  written  as 


v  2  ..-1/2  , .  *  .  .  * 

u  -  i  - J7J  =  ■„  +  €  K  -  >  *2y»)  + 


(€K) 


(52) 


Retaining  the  largest  terms  in  the  differential  equation  (22)  gives 
Laplace's  equation 

*  # 

4>,  ♦  *  +  <{>->  sjt  *  =  0 

2x  x  2yy 


(53) 


The  wall  boundary  conditions  are 


4>*  *<x*.±l,t)  =  T  K1/2  x;:;f"(x0), 


0  <  x*  < 


(54) 


for  the  interior  walls,  and  i\>~  *(x*,±,l,t)  =  0  for  0  <  x*  <  oo  along  the 
exterior  surfaces;  these  boundary  conditions  are  indicated  in  Fig.  3a. 
The  leading-edge  solution  for  u  must  match  with  the  channel  solution  for 
u;  that  is,  the  expansion  of  the  leading-edge  solution  as  x*  -*  oo  ,  for 


-1  <  y*  <  1 ,  is  required  to  agree  term  by  term  with  the  expansion  of  the 

,  ,  2  -2 

channel  solution  as  x  -*■  x^,  shown  in  Eq.  (46a).  Finally,  as  x  "  +  y"'  -*■  oo 

outside  the  channel,  the  velocity  perturbations  must  approach  zero. 

The  solution  to  the  flow  problem  formulated  above  can  be  found 
with  the  help  of  a  conformal  transformation  from  the  complex  z* -plane, 
where  z*  =  x*  +  iy*,  to  a  complex  4-plane.  The  boundaries  shown  in 
Fig.  3a  are  therefore  transformed  to  the  real  axis  in  Fig.  3b,  where 
the  leading  edges  z*  =  ±.  i  have  been  mapped  onto  the  points  4  =  ±.  1. 

The  required  transformation  is 

z*  =  -  (42  -  1)  -  -  lnt,  +  i  (55) 

TT  IT 


If  the  complex  potential  corresponding  to  <j>*  is  F(z*)  =  *3f(4),  the  complex 
velocity  is 

F'(z*)  =  7* (4)  d4/dz*  =  <|>2*x*  -  i  <>*  *  (56) 


A  simpler  preliminary  boundary -value  problem  is  obtained  by  differen- 

j{C 

tiation  with  respect  to  x  ,  since  the  wall  values  of  <j>_  $  *  are  constant. 

wX  y 

The  solution  to  this  problem  is 


.  .  # 

1 


2x  y 


*  - 


K1/2f"(x0)  _  j 

-  In  (1  -  — 


(57) 


Integrating,  and  adding  a  source  term  at  the  origin  4  -  0  so  that  the 
matching  condition  with  the  channel  solution  can  be  satisfied,  one  finds 


finally 


K'/2f,(V  2  1 

♦i*  - s  *2>  =  — r-2-  a  *  d2  - 1)  (i  -  -L) 


oo  ,  , 

2  /  ^  In  (1  -~)d4> 


K 


1/2  1  1 
~  {—  f,,(x0)  +  m2(t)}~ 

n 


U  -  U  -  X 

00 


(€K) 
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2  it  z 


r  -i 

(58) 


where  m^(t)  is  still  to  be  determined. 

J*C 

As  z'"  oo  outside  the  channel  (i.e.  ,  with  arg  z  fixed, 

0  <  arg  z*  <  2tt),  the  transformation  (55)  becomes  z*  =  t,  /it  +  •  •  •  ,  and 
expansion  of  F' (zv)  gives  velocity  perturbations 

v  2  m2(t) 


(59) 


This  expression  describes  the  velocity  due  to  a  source  of  strength  m^ft) 

at  z  =  0,  and  provides  a  matching  condition  for  velocity  perturbations 

in  the  external  flow,  which  could  be  calculated  if  desired.  For  such  an 

outer  solution  the  proper  length  scale  in  the  flow  direction  would  be  the 

duct  length,  so  that  x  =  0(1),  In  the  differential  equation  the  terms  «j> 

and  (Y+l)<j>  <J>  (Y+l)£u.  <b_  are  then  of  the  same  order  (and  all 

X  XX  1  oo  xx 

other  terms  are  of  higher  order)  when  y  =  Ofe  so  that  an  appro- 

priate  y-coordinate  would  be  y  =  (Ke)  '  y.  For  x  =  0(1)  and  y  =  0(1) 
the  largest  perturbations  on  the  uniform  flow  are  described  by  a  velocity 
potential  which  satisfies  Laplace's  equation  in  x  and  y.  For  matching 
with  the  leading- edge  result,  the  solution  would  contain  one  term  which 
represents  a  source  at  x  =  0,  y  =  0;  the  remaining  part  of  the  solutions 
would  be  determined  by  the  prescribed  unsteady  normal  velocity 
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component  at  the  (exterior)  duct  walls  located,  to  the  outer  length  scale, 
at  y  =  0.  The  details  of  such  outer  solutions  are  not,  however,  of 
special  interest  here. 

At  the  edges  z  =  ±  i,  the  leading-edge  solution  for  the  velocity 
has  the  inverse-square-root  singularity  associated  with  flow  around  a 
sharp  edge.  As  z  -*-±.1, 


u  -  u 


IV 


00 


(eK) 


1/2 


m  (t)  +  f"(x  )/ 7T 

~  ±  e  - rrx-  “P  (-  7  arg  (z*  ”  i))  (&°) 

2R„U*Ti|)1/2  2 


A  positive  or  negative  value  of  m^ft)  +  fn(xQ)/iT  implies  an  outflow  or  an 
inflow,  respectively,  around  the  edges.  An  outflow  around  the  edges 
remains  even  if  the  source  strength  m^t)  is  zero,  since  the  flow  ahead 
of  the  duct  anticipates  the  area  decrease  and  accelerates  slightly  before 
reaching  the  duct.  Thus  a  slight  contraction  of  streamtube  area  occurs 
upstream  of  the  duct,  and  the  small  excess  mass  flow  spills  around  the 
edges. 

The  function  m^ft)  is  found  from  the  matching  of  the  leading-edge 
and  channel  solutions.  As  x  *♦  oo  with  -1  <  y  <  1,  the  transformation 
(56)  implies  that  £  -*■  0,  and  the  expansion  as  £  -*•  0  of  the  solution  (58) 
can  be  shown  to  give 

u  -  u  =  T  {f"(x  )  (x*2  -  y*2  +  4)  -  m  (t)}  +  •  •  •  (6la) 

oo  2  u  3  2 

v  =  -  K^2  x*  y*  f'^x^)  +  •  •  •  (6lb) 

Comparison  of  Eqs.  (61a)  and  (46a)  shows  that 
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by  the  imposed  pressure  downstream  of  the  duct,  as  shown  in  Eq.  (44) 
for  cases  in  which  no  special  trailing-edge  solution  is  needed.  For 
example,  the  duct  might  have  zero  slope  and  curvature  at  trailing  edges 
which  are  joined  to  fixed  parallel  walls  downstream;  or  the  flow  might 
exhaust  into  a  plenum  at  x  =  x^,  again  with  f^'frj)  =  f^'(xj)  =  0  and  with 
the  flow  leaving  the  edges  smoothly,  i.  e.  ,  with  the  wall  streamlines 

having  continuous  curvature.  Of  primary  interest  here,  however,  is  the  i 

case  of  choked  flow,  for  which  c^(t)  is  determined  by  conditions  at  the 
throat,  as  anticipated  in  Eq.  (45)  and  as  discussed  below. 

Flow  Details  Near  the  Throat 

At  a  throat,  the  flow  velocity  may  have  a  maximum  (for  subsonic 
flow)  or  a  minimum  (for  supersonic  flow),  or  the  flow  may  accelerate 
or  decelerate  through  the  sound  speed.  In  any  case,  if  |u  -  1  j  «  e 
near  a  throat  at  x  =  0,  then  4>^x  **■  0  as  x  -*  0  and,  from  Eq.  (40a), 
h^  —  oo  as  x  -*  0,  Since  4^  therefore  does  not  remain  small  in  com¬ 
parison  with  <t>jx  near  x  =  0,  it  should  be  anticipated  that  the  channel- 
flow  formulation  may  not  always  be  correct  near  a  throat,  and  that  for 
some  combinations  of  the  parameters  a  local  "inner"  solution  in  a  thin  j 

throat  region  might  be  needed.  In  order  to  understand  what  modifications  I 
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might  be  required,  one  might  consider  a  purely  av  3onic  flow  with  the 
back  pressure  gradually  decreasing  (or  the  undisturbed-flow  Mach  num¬ 
ber  increasing)  until  sonic  speed  is  reached  in  the  throat  region.  This 

approach  is  taken  in  the  following  paragraphs,  and  two  special  cases 

3/2  2 

are  seen  to  arise,  corresponding  to  u  -  1  =  0(€  )  and  u  -  1  =  0(€  ) 

near  the  throat.  In  the  latter  case,  acceleration  to  supersonic  speed  is 
considered  in  detail. 

According  to  the  channel  solution  (41a),  the  (subsonic)  velocity 
near  a  throat  located  at  x  =  0  is 


1  -e  (^J)  '  |Cl  -  f(0)  -|x2f"(0)+  ~)1/2  +  e2j^f"(0)  (}-y2) 

+  i[fi,(0)-fu(0)ly  +  (i-l)  vti  [ci  -f(0)] 


G(0,t)  +  c2(t) 

[2(y+l)]1/2  [Cj  -  f(0) -ix2  f"(0)  +  ...]1/2 

+  • • •  (63) 

2 

where  c.  =  ("Y  +  1)  u.  /2,  from  Eq.  (48).  If  the  flow  is  choked,  then 
1  1 00 

2 

c.  =  f(0);  i.  e.  ,  ("Y  +  1)  u.  /2  =  f (0).  The  second  term  in  Eq.  (63)  is  then 

X  1 00 

1  /  2 

no  longer  small  in  comparison  with  the  first  term  when  x  =  0(e  )  and 

3/2 

u  -  1  =  0(€  ).  Ifa  special  local  solution  corresponding  to  these 

orders  is  assumed  necessary,  one  finds  that  the  differential  equation 
and  boundary  conditions  lead  to  a  problem  formulation  for  x  =  Ofe^2) 


which  is  similar  in  form  to  that  which  leads  to  the  solution  (37a)  for 


<)>jx.  The  local  solution  is  found  to  give 

u  =  1  ■ e3/z  (y+i) 7  {- It*  fU(0)  +  c2(t)  +  G(0,t)}1/2  +  ■"  (64) 

1/2 

where  the  functions  of  t  are  chosen  for  matching  as  x/e  -»  -oo  with 
the  expansion  (63)  of  the  channel  solution  as  x  -*•  0.  It  can  be  seen  that 
c^(t)  +  G(0,t)  >  0  for  subsonic  flow  at  the  throat.  The  result  (64) 
remains  correct  at  x  =  0  and  simply  describes  a  change  in  fluid  accel¬ 
eration  from  positive  to  negative  values  within  a  short  distance 
1/2 

x  =  0(c  )  of  the  throat.  A  convenient  composite  solution,  suggested 

in  Ref.  [17]  for  steady  flow,  is  obtained  by  replacing  the  integration 
constant  c^  with  a  suitable  series  expansion  in  e  : 

2  1/2 
u  =  1  -  €  -  f(x)  +  e  c2(t)  +  e  G(x,t)  +  •  •  •  }  ' 

+  €2{-|f"(x)  (y  -y2)  +  ~  [f|'(x)  -  r'(x)]  y  +[j  -  j)  [cj  -  f(x)]  } 

+  •  •  •  (65) 

The  solution  (41a)  is  recovered  if  the  square  root  is  expanded  for 
|e  c^t)  +  €  G(x,t)|  <<  Cj  -  f(x),  and  if  Cj  =  f(0),  the  local  solution  (64) 
is  obtained  when  x  =  0(e^2  ). 

If  now  Cj  -  f(0)  =  0  and  also  c^(t)  +  G(0,t)  =  0,  either  the  channel 

solution  or  the  composite  solution  implies  that  no  difficulty  now  arises 
1/2 

when  x  =  0(€  ),  but  that  when  x  =  0(e )  the  largest  term  in  u  -  1  is 

2 

0(€  )  and  is  a  function  of  y  as  well  as  x/€ .  This  conclusion  also  follows 
from  consideration  of  the  differential  equation  and  boundary  conditions. 
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The  wall  boundary  condition  (31)  shows  that  ^  =  0(e  x)  as  x  —  0.  In 
the  potential  equation  (22),  it  is  expected  that  (Y  +  1)  4>  <|>  is  no  longer 
small  in  comparison  with  when  the  local  Mach  number  is  sufficiently 
close  to  one.  Since  also  y  =  0(1),  it  follows  from  equating  the  orders  of 
these  terms  that  a  local  solution  satisfying  a  nonlinear  differential  equa- 
tion  near  the  throat  may  be  needed  when  x  =  0(e).  Coordinates  x  '  and 
y*  are  defined  by 


x*  = 


(y  +  i)1/2  [-f"(0)31/2e 


* 

y  =  y 


(66) 


and  the  perturbation  potential  is  represented  in  the  form 


4>  =  (Y+l)1/2  [-f"(0)]3/2  €3  4>*(x*,y*,t)  + 


(67) 


The  largest  terms  in  the  differential  equation  and  boundary  conditions 
give 

^Vy*  =  *2**  *2fc*x*  (68a) 


* 

+2y* 


.  f  i'(0)  -  f  "(0)  . 

(x*.  ii.t)  =  |±i  -  — TnoT~ } x* 


{ 


+  <  T  1  + 


Glx(°»t)  -  Gux(0>t) 
2  G  (0,t) 


K 


(68b) 


where 


x*  =  -  (Y  +  l)"1/2[-fM(0)]"3/4  Gx(C, 


-3/2 


(69) 


Of  particular  interest  is  the  solution  which  describes  acceleration 
of  the  flow  from  subsonic  to  supersonic  speeds.  As  for  steady  flow,  4>-j! 
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is  found  as  a  polynomial  in  x*  and  y*,  but  now  some  time-dependent 
terms  appear: 


#  1  .  jfc.2  1  ,  S{S  ,  #2  1  .  1  .  1  #4  , 

*2s2tx  ~xl}  +l(x  "Xl)(y  -  3)  +  TI  (2y  _y  > 


1  7  2 


3  7  12  '2 


-  “{f;'(0)  -  f»(0)Hx*y*  (-|y*3  -  y*)} / f"(0) 

"2  {G,x(0*t,+  0^(0.*)  >  x*y*/Gx(0,t) 


In  the  original  variables  the  velocity  components  for  x  =  O(e)  are  then 

.  2  (x[-f»(0)]1/2  Gx(0.t)  i  .  2  1 

U  =  1  +  €  { - —  + - Tn - 1777  ■  7  f  (0)  (y  ~  7* 

l(v+i)1/2€  (y  +  1)1/2  [-f"(0)]1/2  2 


\  [i"(0)  -  r*(0)]  y| 


(71a) 


v  =  e3  { -  f"(0)  J  y  +  Gx(0,t)y  +  ^(Y+l)1/2  [-f"(0)]3/2  (y3  -  y) 
+  7  [f ;’(0)  -  f"(0)]  [7  +  7  (Y  +  1)1/2  [-f"(0)]1/2  (y2-  1)] 

b  *  VI  6  b 


+  7  [G.  (0,t)  +  G  (0,t)3  }  +  ••• 

2  lx  ux 


(71b) 


The  sonic  line,  defined  by  u  =  1,  is  located  at 

+^,v+ 1,U2  t'£"'°"l/2(^/) 


.  ..i/2  f;'<o>-£>  ) 

(Y  +  1 )  J77  Y  >  + 

2  [-f"(0)]1/2  ) 


The  second  and  third  terms  give  the  steady-state  result,  and  the  first 
term  is  seen  to  be  the  value  of  x  at  which  the  instantaneous  channel 
width  is  a  minimum.  That  is,  the  sonic  line  is  found  to  have  the  appro¬ 
priate  quasi- steady  location  at  each  instant.  This  conclusion  follows 
only  because  of  the  particular  choices  made  here  for  the  orders  of 
magnitude  of  t  and  a.  For  the  other  problems  discussed  in  subsequent 
sections,  the  sonic  line  does  not  remain  close  to  a  quasi-steady  location. 

The  expansion  (63)  of  the  channel  solution  as  x  -*  0  can  now  be 
matched  with  the  throat  solution  (71a)  to  give  the  results  anticipated 
earlier  in  Eqs.  (38)  and  (45): 

Cj  =  f(0)  ,  c2(t)  =  -  G(O.t)  (73a, b) 

Also  c  =-7(V  +  l)u2  ,  from  Eq.  (48),  and  so  u  =  -  2^  2  f(0 )/  (V  +  l)1^. 

1  2  loo  loo 

That  is,  at  the  minimum  Mach  number  for  choking 

u  =  1  -  e  21/2  f(0)/(y  +  1)1/2  +  0(e2)  (74) 

00 

The  function  c^t)  is  found  here  from  consideration  of  the  flow  details 
near  x  =  0,  and  the  result  of  course  agrees  with  Eq.  (45),  which  was 
obtained  from  the  requirement  that  <t>2x  remain  finite  as  x  -*  0.  Com¬ 
parison  of  the  channel  and  throat  solutions  (41)  and  (71)  shows  that, 

for  acceleration  from  subsonic  to  supersonic  speeds,  the  channel  solu- 

2 

tions  for  u  and  v  do  in  fact  remain  correct  to  order  c  in  the  thin 
region  x  =  O(e)  which  includes  the  throat.  A  special  local  solution 
would  be  needed  here  only  for  flows  which  remain  either  subsonic  or 
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supersonic  when  x  =  O(e),  but  which  approach  so  close  to  the  sound 

1  i  2 

speed  that  the  minimum  (nonzero)  value  of  |u  -  1  |  is  0(e  ).  In  the 
present  case,  the  ratio  {c2(t)  +  G(x,t)}/4>^x  is,  however,  indeterminate 
at  x  =  0  and  must  of  course  be  evaluated  properly  there.  The  combina¬ 
tion  Cj  +  e  c^  appearing  in  the  composite  solution  (64)  is  equal  to 
f(0)  -  e  G(0,t),  which  is  just  the  (relative)  area  change  at  x  =  0.  This 

earlier  composite  solution  can  be  modified  to  include  the  case 

2 

c^(t)  =  -  G(0,t)  by  addition  of  a  suitable  0(c  )  term  inside  the  square 
root: 

2  2  ,  1/2 

2  6  Gx  (0’t)) 

u  =  1  +  e  (— j  ^  -f(x)  +  €c2(t)  +  €G(x,t)+  — -[7FTdjy f 

+  €2{|f"(x)  (j  -  y2)  +~[f;'(x)  -  r*(x)]y+  (7-7)777  [cJ  -f(x)]} 

+  ...  (75) 

M  2 

This  new  composite  solution  is  correct  to  order  e  both  near  and  away 

from  the  throat,  since  it  may  be  expanded  for  x  =  O(e)  to  give  the  throat 

solution  and  for  x  =  0(1)  to  give  the  channel  solution.  Moreover,  the 

case  c^'t)  +  G(0,t)  >  0  is  also  included,  since  then  Eq.  (64)  is  recovered 

1/2 

by  expansion  of  Eq.  (75)  for  x  =  0(e  ).  Thus  Eq.  (75)  is  sufficient  for 

purely  subsonic  flow,  with  c ^ (t )  +  G(0,t)  >  0,  or  for  acceleration  from 

subsonic  to  supersonic  speed,  with  c^t)  +  G(0,t)  =  0,  and  special  con- 

1/2 

siderations  are  no  longer  necessary  for  x  =  0(e  )  or  x  =  0(e). 


Velocity  and  Position  of  the  Shock  Wave 

Just  as  in  the  case  where  the  back  pressure  alone  oscillates 

[9-12]  ,  the  position  of  the  shock  wave  at  any  instant  is  found  by  deriving 

first  the  expression  for  the  instantaneous  velocity  of  the  shock  wave  and 

integrating  this  equation.  The  velocity  is  found  by  satisfying  the  jump 

conditions  across  the  shock  wave.  It  was  shown  earlier  that  the  flow  is 

steady  in  first  order.  If  the  velocity  of  the  shock  wave  were  O(e),  then 

unsteadiness  would  be  induced  in  the  O(e)  velocity  terms  downstream  of 

2 

the  shock  wave,  so  it  is  clear  that  for  the  conditions  chosen  u  =  0(e  ); 

s 

this  is  consistent  with  the  order  of  Ug  found  using  Eqs.  (7b)  and  (32b). 

If  the  expansion  for  u  given  in  Eq.  (5a)  and  the  value 
H  =  OY  +  l)/2(y  -  1)  are  substituted  into  Eqs.  (30c, d),  it  is  seen  that  the 
jump  conditions  across  the  shock  wave  are 


(76a) 


4  Us 

U2U  '  Uld  Ulu  +  (Y  +  1)  (£2  ] 


(76b) 


It  is  clear  from  Eqs.  (41a)  that  if  the  solution  for  u^  with  the  positive 
sign  (supersonic  flow)  is  associated  with  the  flow  upstream  of  the  shock 
wave  and  that  with  the  negative  sign  (subsonic  flow)  is  associated  with 
the  flow  downstream  of  the  shock  wave,  the  jump  conditions  on  Uj  are 
satisfied.  On  the  other  hand,  because  of  the  dependence  of  u^  on  y,  it 
is  clear  also  that  the  second-order  jump  conditions  cannot  be  satisfied 
by  the  channel-flow  solutions;  additional  calculations  are  required. 
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Because  the  second-order  jump  conditions  are  not  satisfied,  it  is 

necessary  to  consider  a  thin  inner  region  behind  the  shock  wave,  within 

which  some  adjustment  of  the  flow  must  take  place.  The  solutions  in 

1/2 

this  inner  region,  which  turns  out  to  have  an  extent  Ax  =  0(e  ),  must 

satisfy  the  jump  conditions  across  the  shock  wave  and  the  boundary  con¬ 
ditions  at  the  walls,  and  must  match  with  the  outer  channel  flow  solu¬ 
tions  in  the  appropriate  limit.  The  problem  here  is  complicated  by  the 
fact  that  the  shock  wave,  and  therefore  the  inner  region,  is  moving  such 
that  Axg  =  0(1).  Solutions  for  such  an  inner  region  enclosing  the  shock 
wave  were  first  found  for  steady  flow  [  1 7] .  This  solution  was  extended 

[9]  to  cover  unsteady  flow  with  small-amplitude  motion  of  the  shock 

-1  2 

wave,  with  Ax  =  Ofc)  and  r  =  0(€  ),  so  that  u  =  0(e  ).  In  this  analy- 

s  s 

sis  second-order  oscillations  in  the  back  pressure  caused  the  unsteadi¬ 
ness  in  the  flow,  the  walls  being  stationary.  Next,  it  was  shown  in 

Refs.  [15]  and  [  1 2]  that  for  impressed  pressure  fluctuations  of  the  same 

-2  2 
order,  but  with  t  =  0(e  )  so  that  Ax  =  0(1)  and  u  =  0(e  ),  the  same 

s  s 

form  of  solution  resulted  in  the  inner  region.  Because  no  time  deriva¬ 
tives  remain  in  the  governing  equations  for  this  region,  to  the  order 
considered,  the  differences  in  the  solution  arise  only  from  the  differ¬ 
ences  in  the  matching  conditions  to  be  met;  for  t  =  0(e“*)  the  shock  is 
stationary  and  the  incoming  velocities  evaluated  at  the  shock  are  con¬ 
stant  to  lowest  order  while  for  t  =  0(e  ^),  the  position  of  the  wave  and 
the  incoming  velocities  are  functions  of  time  to  lowest  order.  As  a 
result,  the  equations  for  the  velocity  of  the  shock  wave  have  the  same 


dependence  upon  h^,  etc.  in  the  two  cases,  but  in  each  case  these 

2 

functions  depend  differently  upon  time.  Also  u  =  e  dx  /dt  for 

S  SI 

-1  2  -2 
t  =  0(e  ),  whereas  u  =  e  dx  rt/dt  for  t  =  0(e  ),  as  seen  from  Eq. 

s  su 

(7b). 

The  analyses  referred  to  in  Refs.  [l7,  9,  15,  12]  all  involved 
symmetric  channels.  However,  Richey  [18]  considered  the  case  of 
asymmetric  walls  also,  for  t  =  0(e  *),  and  showed  that  a  simple  trans¬ 
formation  of  the  y  variable  in  the  inner  region  could  be  used  to  relate 
the  solutions  for  symmetric  and  asymmetric  walls.  His  outer  channel- 
flow  solutions  have  the  same  form  as  those  found  here  to  second  order, 
even  though  t  is  of  different  order  in  each  problem,  because  in  both 
cases  the  flow  is  steady  in  first  order  and  the  walls  are  stationary  to 
second  order.  In  view  of  the  above  remarks,  it  is  clear  that  in  the 
inner  region  behind  the  shock  wave,  those  solutions  found  before,  modi¬ 
fied  for  asymmetric  wall  shapes  following  Richey  [l 8 ] ,  hold  in  the 
present  case.  The  composite  solutions  for  u  and  v,  formed  from  the 
inner  and  channel-flow  solutions,  and  the  equation  for  ug  are  as  follows: 

u  = 1  ±  (ci-f)  +  e2[-f"  V +  +  h2x +  ^2xJ  + (77a) 

v  =  e2  [  -f  y  +  g2]  +  €  5//2  t2y  +  *  '  *  (77b) 

£  =  0  for  x  <  x  ,.  +  •••  (77c) 

su 


40 


A  32  V(V+1)U1U 
^  = - 3 - 

IT 


nt  i 

|g20  f  ^  exp(-M5/V(V+l)aIo)Sinaf- 

'  n=l,3,5,*  •  •  n 


-  V  s  ^ 

n=2,4, 6,  •  •  •  n 


exp  (-  ^ /V^lhT^)  cos 


for  x  >  x  + 
sO 


(77d) 


x  =  (x  -  xg0)  g"1/2  (77e) 

4  Us  4k  ^sO  ,,  .  ..  .  f0  2 

(Y  +  1)  e2  =  (Y+  1)  dt  "  (h2x)d  h2x  u  "  3  “  Ulu  (77f) 


where  Uj^  is  given  by  Eq.  (37a),  with  the  positive  sign,  evaluated  at 

x  =  x  .  +  • .  •  and  is  thus  a  function  of  time,  and  where  the  subscript  0 
sO  r 

indicates  that  the  function  is  evaluated  at  x  =  xsq-  The  equation  for 

2 

ug/€  shows  only  the  first  term  in  the  expansion  since  that  is  all  that 

is  desired  here.  The  composite  solutions  for  u  and  v  are  uniformly 
2  5/2 

valid  through  0(e  )  in  u  and  0(€  )  in  v  throughout  a  channel  in  which 

a  shock  wave  occurs,  as  long  as  the  throat  or  leading-edge  regions  are 
avoided. 

If  Eqs.  (40a)  and  (76a)  are  used  to  find  (Vi2x)  and  (h^,,)  and 

d  u 

the  results  are  substituted  into  Eq.  (77f)  ,  one  finds,  finally,  the 
differential  equation  for  the  position  of  the  shock  wave  to  lowest  order, 


TTn  ~dT  =  ’  7"~rn -  ^  (v  + 1)ui3  '  c2  +  czJ  (78) 

(y  +  l)  dt  (y+l)Uju  3  lu  2u  2d 


where  c.  and  c_  ,  are  the  values  of  c_(t)  used  upstream  and  downstream 
2u  2d  2 
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of  the  wave,  respectively.  If  the  example  discussed  previously  is 


considered,  where  the  flow  is  subsonic  entering  the  duct  and  acceler¬ 
ates  to  supersonic  velocity  upstream  of  the  shock  wave,  with  pressure 

oscillations  impressed  upon  the  flow  at  x  >  x  ,  then  Eqs.  (44)  and 

X  s 

(45)  hold  for  cOJ  and  c  ,  respectively;  then  Eq.  (78)  becomes 
£.d  lm 


4k  dXsO 


(V  +  1 )  dt 


77 7TT -  [~(V  +  l)u.3  +  N(t)] 

(y  + 1)  u.  3  lu 

lu 


(79a) 


N(t)  =  C2d  }  +  °d(t)  +  G(0,t)  "  G(xi*t) 


(79b) 


Although  analytical  solutions  may  be  found  for  special  wall  shapes  and 
forcing  functions,  the  calculation  of  generally  will  involve  numeri¬ 
cal  integration  of  this  equation. 


Motion  of  the  Shock  Wave  Upstream  of  the  Throat 

According  to  Eq.  (78)  or  (7  9a)  the  shock- wave  velocity  dx  /dt  —  oo 

S  v 

if  the  upstream  fluid  velocity  u^  -*•  0,  as  occurs  if  the  shock  wave 
approaches  the  throat  as  the  time  t  approaches  a  value  t,,..  Expanding 
Eq.  (79a)  as  xg^  -*■  0  and  t  —  t*  ,  and  then  integrating. 


4k  dXs0 
V  +  1  dt 


N(tJ 

(y  +  D1/2  [-f"(0)]1/2xs0 


fY  +  1)1/2  N(t*) 

- 77 j  (t„.  - 1)  +  ... 

2k  [-f'MO)]1^ 


(80a) 


(80b) 


where  N(t;J  = 


C2d^*^  +  G^,t:!-')'  *n  accordance  with  the  definition  (79b). 
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If  an  initial  shock-wave  position  has  been  specified,  the  complete  solu¬ 
tion  for  x  rt(t)  can  be  regarded  as  known,  and  therefore  t .  is  known, 
su  "■ 

When  the  shock  wave  is  near  the  throat,  a  different  form  of  solution 
is  required.  It  is  shown  below  that  the  shock-wave  velocity  increases 
in  order  of  magnitude  as  it  passes  through  the  throat  region  and 
becomes  still  larger  as  it  moves  further  upstream  in  the  duct.  The 
end  results  obtained  from  the  analyses  of  this  and  the  preceding  sec¬ 
tion  are  expressions  for  the  shock-wave  velocity  in  terms  of  position 
at  all  points  of  an  accelerating  channel  flow  which  is  subjected  to  pre¬ 
scribed  small  oscillations  in  wall  shape  or  back  pressure. 

The  composite  solution  (65)  written  for  the  velocity  in  the  region 
downstream  of  the  shock  wave  gives 


2  1/2  1/2  2 
u=i-€(— )  {f(0)  -  f(x)  +  ec2d(t)  +  e  G(x,t)}  '  +  Ofe  )  (81) 


where  in  general  c2d(t)  /  c^t)  and  so  c2d^  N(t)  +  0. 

As  the  shock  wave  moves  upstream,  the  perturbation  in  fluid  velocity 

3/2 

behind  the  shock  wave  found  from  this  solution  is  0(c  )  when 

x  =0(e1/2): 
s 


ud  =  1  -  e3/2  (^~)  {-  f-  f  "(0)  +  N(t)}1/2  +  0(e2) 


(82) 


The  velocity  just  upstream  of  the  shock  wave,  since  c.  (t)  +  G(0,t)  =  0 
and  f"(0)  >  0,  is 


,/p  p  l/Z 

.3/2  ,  2 _  ,  ,  x 


uu=lt£"‘l7ri>  '“Vi1 


(83) 


43 


The  shock-wave  jump  conditions  (30d,  e),  with  H  =  —  (Y  +  1)/  (V  -  1),  give 


^  {(«u  -  1)+  (ud  -  1)}  + 


The  shock-wave  speed  u  *  dx  /dT  is  therefore  0(6  ),  and  the 

s  s 
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shock  wave  remains  within  a  distance  x  =  0(e  )  of  the  throat  for  a 

s 

- 1  2 

time  interval  AT  =  0(€  ).  Since  t  =  ke  T,  the  interval  in  t  is  small, 

namely  At  =  0(e).  The  shock  wave  therefore  passes  through  the  thin 

1/2 

region  x  =  0(€  )  at  essentially  a  constant  value  of  t.  For  matching 

s 

with  the  downstream  solution,  as  shown  below,  this  constant  has  the 

value  t^  obtained  in  E q.  (80b),  Hence  t  is  replaced  by  t*  +  O(e)  in  the 

function  N(t),  with  u  expected  to  be  obtained  as  a  function  of  (t  -  t.,.)/e, 

S  ' 

2 

Since  u  =  k€  dx  /dt,  Eq.  (84)  becomes 
s  s 


Z*1,  3/2  (V+l)i/2  1  x 

ke  IT  =  e  — 3 — {7 


H72 


[-  -j-  f"(0)  +  2  N(t*)]  1/2  |  + 


The  shock-wave  position  therefore  has  the  form 


xa  =  €l,c  [-2N(t*)/f"(o)]1/2  eg;j(t:::>  + 


(86a) 


t*  =  4£  (V+D1/2  [-f"(0)]I/2  (t  -  t.,.)/e 


(86b) 


Integration  gives  an  implicit  solution  for  £ g  j  (t  ): 


(87) 


2t*  =  !-(!  +  4*f)1/2  {g*,  +  (1  +  g*2)1/2> 


-  In  {g*  +  (1  +  g*?)1/z} 
si  si 

As  gg  j  •*  oo  and  t  -*■  -oo  ,  this  solution  gives 

I*  =  (-t*)1/2  +  ...  (81 

S  1 

The  matching  as  x  .  *♦  0  and  g*  -*•  oo  requires  that  the  largest  terms 
su  si 

of  Eqs.  (80b)  and  (88)  be  identical,  when  written  in  common  notation. 

It  is  this  condition,  anticipated  above,  which  requires  that  the  time  in 
Eq.  (86b)  be  measured  relative  to  the  value  t#  obtained  from  the  down¬ 
stream  solution. 

Next,  as  the  shock  wave  moves  upstream  from  the  throat,  for 

g*  -*  -oo  and  t*  -*  oo  ,  it  is  found  from  the  solution  (87)  that  In  (-2g*  ) 
si  si 

=  2  (t*  -  "T )  +  •  •  *  ,  or 
4 


]/2 

N(t*)  I1/2  2  '  k€ 

-2ef”(0)J  * 


In  (l/€ ; 


(V+l(/2[-£"(0)]I/J 


When  x  =  0(1),  the  channel  solution  (3'a)  has  the  subsonic  value 
s 

(minus  sign)  both  upstream  and  downstream  of  the  shock  wave,  so  that 

Uj^j  =  and  the  shock-wave  velocity  from  Eq.  (84)  becomes 

u  =  —  0Y+  1)€  u,  +  •  •  •  ,  That  is,  using  also  Eq.  (21b),  to  order  e 
s  Z  lu 

the  shock  wave  travels  upstream  at  a  speed  u  =  u  -  a  +  .  . .  ,  and  the 

order  of  magnitude  of  u  has  now  increased  further,  so  that  u  =  0(€). 

s  s 

The  solution  for  the  position  of  the  shock  wave  as  it  moves  upstream 
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from  the  throat,  with  x  =  0(1),  must  have  the  form  (89)  as  x  -*  0. 

s  s 

That  is,  the  solution  for  x  =0(1)  must  match  with  the  solution  for 

s 

1/2 

x  =  0(e  )  at  earlier  times.  This  condition  suggests  a  representation 

s 

of  the  shock-wave  position  in  terms  of  still  another  time  coordinate,  as 
fo  llows : 

Xs  =  Xs0  +  “  '  (90a) 


t*  _  *  "**  _  In  (1/e) _ 

k€  (Y+l)1/2  [-f"(0)]1/2 


(90b) 


Substitution  in  the  expression  (84)  for  the  shock-wave  velocity  gives 


dx  + 
_ sc 

dt+ 


Y+1I/2  1/2 

=  -  l1^)  { f (0 )  -  f(x+0)}1/2 


(91) 


Since  f(0)  -  f(x+  )  ~  x+2  f"(0)  as  x+_  -*  0,  the  solution  as  x+.  ■*  0 

sO  2  sO  sO  su 

and  t  *♦  -oo  has  the  form 


\  (Y  +  l)1/2  [-f"(0)]1/2  t+ 


r  N(t„)  -iI/2 

|_  -  2  e  fM(0)J 


+ 


(92a) 

(92b) 


where  A  is  an  integration  constant.  Comparison  with  Eq.  (89)  shows 
the  exponents  are  identical,  and  the  requirement  that  the  solutions  must 
match  has  been  used  to  give  the  value  of  A.  This  result  provides  an 
initial  condition  for  the  differential  equation  (91). 

The  simplest  special  case  arises  for  a  duct  of  length  2x^,  with 
Xq  =  -  Xj,  having  symmetric  parabolic  walls  described  by 
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i 


(93) 


f(0)  -  f(x)  =  f(0)  — 
xi 


where  f(£x.{  =  0  and  the  choking  condition  (74)  gives  f (0 )  =  -y  (V  +  I)u  ^  . 

Integration  of  the  shock-wave  velocity  (91)  shows  that  the  solution  (92), 

with  f"(0)  =  -  (V  +  1)  u,^  /xf,  remains  correct  for  the  entire  shock-wave 

1  00  1 

motion  from  the  throat  x  =  0  up  to  the  entrance  x  =  -  x  .  The  time 
required  for  passage  of  the  shock  wave  through  this  region  is  found  to 


k  x.  e  2  (Y  +  1 )  u.  e 

.  _ 1 _  .  _ loo 

*  "  (Y  +  D  (-uloo)  m  £N(g 


The  shape  (93)  implies  wedge-shaped  leading  and  trailing  edges,  since 
f1  (i.Xj)  ^  0.  A  simple  shape  having  cusped  edges  is  given  by 


2  irx 


f(0)  -  f(x)  =  f(0)  sin 


1  2 

where  again  f(i  x  )  =  0  and  f(0)  =—  (Y  +  1 )  u.  .  Eq.  (91 )  can  be  inte- 

i  L.  1  00 

grated,  and  the  integration  constant  is  found  by  matching  with  Eq.  (92) 

1  Z  2  2 

as  x*  •*  0,  now  with  f”(0)  =  -  —  it  (Y  +  1)  u,  /x,  .  The  shock -wave 
sO  4  loo  1 

position  and  the  time  required  for  the  shock  wave  to  move  upstream 
through  the  front  half  of  the  duct  are  found  as 


_  N(V 

4xl  8(Y+l)eu.2 

L  loo  J 


(Y  +  1  )  IT  (-U.  )t 
1/2  _ 1°° 

4x, 


2xj  k€ 

4  "  **  =  (Y  +  1 ) it  (-u1qo)  ln 


8(Y+  Deuj 
£  N(tJ 


(96a) 


(96b) 


IWfe'LW  J  ^9*  < 


The  results  (94)  and  (96b)  for  different  shapes  show  identical  dependence 
on  the  parameters;  only  the  numerical  factors  are  different. 


While  the  dependence  of  the  solutions  on  the  parameters  is  shown 
explicitly  in  these  special  cases,  applications  for  arbitrary  shape  and 
motion  of  the  walls,  and  arbitrary  oscillations  in  back  pressure,  will 
require  numerical  solution.  In  the  preceding  derivations  the  shock- 
wave  velocity  u  is  found  from  Eq.  (79a)  for  x  =  0(1)  downstream  of 

S  8 

1/2 

the  throat,  from  Eq.  (85)  for  xg  =  0(e  )  near  the  throat,  and  from 

Eq.  (91)  for  xg  =  0(1)  upstream  of  the  throat.  For  numerical  calcula¬ 
tions  some  procedure  is  needed  to  indicate  what  solution  is  to  be  used 
at  each  particular  location.  It  is  probably  most  convenient  to  replace 

these  three  representations  by  two  composite  expressions,  one  for 

2 

x  i  0  and  one  for  x  SO,  and  both  uniformly  valid  to  order  e  .  This 

3  3 

formulation  can  be  accomplished  for  x  <  0  by  adding  the  downstream 

and  throat  representations,  and  then  subtracting  the  common  part  found 

1/2 

for  x  «  1  and  x  »  €  ;  for  x  <  0  the  upstream  and  throat  solutions 

s  s  s 

are  added,  and  the  common  part  for  |x  |  «  1  and  |x  |  is  sub- 

s  s 

tracted.  The  result  for  x  >  0  is 

s 


u  = 
s 


2  y  *  v+i  1/2  f  21/2 

-  e  T  [f (0 )  -  f(x  )]  +  V  N(t)  { - r 

\  x  [-f”(0)] 
s 


4  '  2 


/2 


[f (0)  -  f(x_)f/< 


Jte3/2a±^jjls_  t-r-wU1'* 


-  [  -~  f"(0)  +  2  N(t)] 


1/2  j 


(97a) 
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and  for  x  <  0 
s 

u  =-  €  (-1)  '  [f(0)-f(x  )]1/2  -  e3/2  (-— j  ■■  [-  £"(0)  +  2  N(t)]1 7 2 

Sc.  S  4  t 

-€7{V+1)1/2x  [-f"(0)]1/2+  ...  (97b) 

4  s 


where  u  =  dx  /dT  =  k  €  dx  / dt. 
s  s  s 


Discussion 

The  solutions  presented  here  hold  for  unsteadiness  caused  by  any 
combination  of  impressed  oscillations  of  the  wall  shapes  or  the  back 
pressure.  As  a  result,  they  allow  easy  comparison  of  the  effects  of 
one  as  opposed  to  the  other.  Although  the  results  presented  here  for 
the  channel  flow  have  a  superficial  resemblance  to  those  given  pre¬ 
viously  [15]  for  a  study  of  channel  flow  with  only  oscillations  in  back 
pressure,  there  are  some  significant  differences,  as  will  be  see...  In 
addition,  the  earlier  analyses  have  been  extended  here  to  provide 
detailed  solutions  for  the  leading -edge  and  throat  regions  and  for  the 
shock-wave  motion  upstream  of  the  throat. 

2 

The  velocity  of  the  shock  wave,  u  =  k  e  dx  „/dt  +  .  .  .  ,  and  its 

s  sO 

position  XgQ  are  given  by  Eq.  (79a)  and  its  integral,  respectively.  It 

is  seen  from  Eq.  (79b)  that  the  forcing  function  for  the  shock  motion 

N(t)  could  have  the  same  form  for  either  oscillations  in  back  pressure 

or  in  wall  shape;  i.  e.  ,  either  G(x,t)  =  G(0,t)  =  G(x1ft)  =  0  with  G 

l  d 
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equal  to  some  function  of  time,  or  =  0  with  G(x,t)  such  that 
G(0,t)  -  G(Xj,t)  equals  the  same  function  of  time.  Hence,  insofar  as 
the  motion  of  the  shock  wave  is  concerned,  one  can  find  an  equivalence 
between  wall  and  back-pressure  oscillations.  This  does  not  apply  to 
the  velocity  and  pressure  distributions,  however,  as  shown  by  the 
equation  for  h.,  (Eq.  (40a)),  which  appears  in  u^  and  P^.  No  matter 
what  the  boundary  conditions  may  be,  G(x,t)  depends  upon  x  as  long  as 
a  pitching  motion  of  the  walls  is  involved,  so  that  c^(t)  +  G(x,t)  also 
depends  upon  x;  on  the  other  hand,  if  only  oscillations  in  back  pressure 
occur,  G  =  C  and  c.,(t)  depends  only  upon  G  (t)  and  is  independent  of  x. 
This  point  will  be  illustrated  later  in  an  example  in  which  the  differences 
are  calculated  for  the  flow  upstream  and  downstream  of  the  shock  wave. 

It  is  of  interest  to  note  that  neither  the  solutions  for  the  motion  of 
the  shock  wave  nor  those  for  the  distributions  of  velocity  and  pressure 
depend  upon  x^,  the  center  of  rotation  of  pitching  motion.  This  is 
shown  in  Eq.  (79b)  for  the  shock  motion;  the  same  result  can  be  demon¬ 
strated  for  the  velocity  and  pressure  distributions  by  forming  c^  +  G 
upstream  and  downstream  of  the  shock  wave,  using  Eqs.  (44)  and  (45) 
for  c  and  c0  .  In  all  cases,  the  differ ence  between  G  functions,  each 
evaluated  at  a  different  point,  is  found;  hence,  the  x^  (3(t)  term  cancels. 
The  net  result  is,  then,  that  because  it  is  evidently  area  differences 
which  are  important  in  determining  flow  conditions,  the  center  of 
rotation  is  not  important. 
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All  of  the  solutions  presented  so  far  have  been  discussed  in  terms 
of  application  to  ducts  in  which  the  entire  upper  and  lower  walls  are 
oscillating.  However,  they  may  be  applied  also  to  the  case  corre¬ 
sponding  to  the  flow  through  a  channel  with  stationary  walls  but  with  a 
separated  boundary  layer  downstream  of  the  shock.  That  is,  we 
imagine  now  an  equivalent  wall  shape  which  is  stationary  (Gf  =  G  =0) 

for  x  <  x  say,  where  x  >  x  ,  but  which  is  hinged  at  x  =  x  and  oscil- 
c  c  s  c 

lates  for  x  >  x^.  It  is  seen  that  the  solutions  for  the  velocity  compo¬ 
nents  and  thus  the  pressure,  temperature,  and  density  present  no 
problems  in  this  application  and  will  not  be  pursued  further  here.  How¬ 
ever,  the  equation  for  the  velocity  of  the  shock  wave  exhibits  interesting 
behavior  which  is  worthy  of  discussion. 

The  equation  which  holds  for  dxgQ/dt  is  Eq,  (78),  where  now  the 

flow  upstream  of  the  shock  wave  is  steady,  but  u,  (x  1  varies  with 

lu  sO 

time  because  it  is  evaluated  at  x  _(t).  If  we  consider  the  case  where 

sU 

the  walls  are  oscillating  as  discussed  above  and  the  back  pressure  is 

oscillating  at  some  point  x  =  x,  ,  now  with  x  >  x  ,  and  where  the  flow 

t  1  c 

is  choked  at  the  throat,  then  this  equation  reduces  to  Eqs.  (79)  with 
G(0,t)  =  0.  Thus , 


4k  ^sO 
(Y  +  1 )  dt 


1  r2l 
(Y+l)ulu  L  3 


(Y  +  1 ) u 


3 

lu 


+  c 


(s) 

2d 


+  Gd(t)  -  G(Xj  ,t)] 


(98) 


This  shows  that  the  effect  of  the  wall  flapping  at  some  point  downstream 
of  the  shock  wave  has  the  same  effect  on  the  position  of  the  shock  wave 


as  an  oscillation  in  back  pressure.  Moreover,  it  is  seen  that  it  is  the 


area  change  at  the  "end"  of  the  duct,  x^,  which  is  important  insofar  as 
the  position  of  the  shock  wave  is  concerned,  since  if  G(x^,t)  =  0,  there 
is  no  effect  whatsoever;  moreover,  the  shape  of  the  oscillating  wall 
between  x  =  x^  and  x  =  x^  evidently  has  no  effect  on  the  position  of  the 
shock  wave  in  any  event.  Finally,  in  this  case,  the  motion  of  the  shock 
wave  is  dependent  upon  x^,  the  center  of  rotation  of  this  flapping  wall. 

Example  Calculations 

It  has  been  shown  that  for  the  orders  of  the  parameters  associated 

with  problem  1,  the  effects  of  the  oscillating  walls  and  back  pressures 

are  interchangeable  insofar  as  velocity  and  position  of  the  shock  wave 

are  concerned,  but  not  in  the  pressure  and  velocity  distributions.  In 

these  distributions,  the  differences  between  the  two  cases  are  seen  in 

the  second-order  terms,  the  first-order  terms  being  independent  of 

time.  As  an  illustration  of  this  feature  of  the  solutions,  reference  is 

made  to  the  results  given  in  Ref.  [  1 5 ]  ;  in  that  analysis,  only  oscillations 

in  back  pressure  were  considered,  the  walls  being  stationary.  For 

example,  Figs.  4a,b,c  and  5  (Figs.  6a,b,c  and  7  in  Ref.  [  1 5 ]  , 

reproduced  here  for  easy  reference)  show  isotachs  and  wall  and  center- 

line  pressure  distributions  at  various  times  and  u  and  x  as  functions 

s  su 

of  time  for  the  following  conditions,  where  the  notation  has  been  changed 
to  agree  with  that  in  the  present  report,  and  where  it  should  be  noted 
that  the  calculations  were  made  for  a  symmetric  channel,  so  that  f  =  ff  =  f 


f(x) 


0  <  x  <  1 


27 

13 


(x-2)4 
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13 


(x-2)3 


1  <  x  <  2 


(99) 


=  0 


x  >  2 


G.(t)  =  3  (V+ 1)  sin  (2t) 
a 


t  >  0 


=  0 


t  <  0 


€  =  0.  1 

T  =  100 


V  =  1.4 


-  -y  (V  +  1 )  (u 


(s) 

lu 


3 

) 


=  0 


(S  )  (g  )  / g  \ 

Here  u.  =  u,  (x  )  where  x  .  is  the  location  of  the  shock  wave  in 
lu  lu  sO  sO 

steady  flow,  for  which  =  0.  Since  no  oscillations  in  wall  shape  were 
considered,  G.  =  G  =  G  =  0.  Thus  the  values  of  c?(t)  to  be  used  in 
h2x  (and  thus  in  u^,,  p^»  etc.  ),  upstream  and  downstream  of  the  shock 
wave,  and  the  equation  for  dxgQ/dt  (Eq.  (79a))  are,  for  this  case. 


C2u  +  G  *X,t^  =  c2^  =  0 


(100a) 


/g  )  (g  \ 

c2d  +  G(x,t)  =  c2d  =  c2d  +  Gd(t)  =  c2d  +  3  (y+1)  sin  (2t)  (100b) 


(Y 
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where  k  =  (t  e  )  .  It  is  important  to  note  that  in  this  case,  as  illustra¬ 

ted  by  Eq.  (100a),  the  flow  upstream  of  the  shock  wave  is  steady. 
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If  we  now  consider  flow  in  a  channel  with  symmetric  oscillations 

of  the  wall,  so  that  G,  =  G  =  G,  and  with  no  oscillations  in  back  pres- 

t  u 

sure,  so  that  G  =  0,  we  see  that  Eqs.  (79)  and  (100c)  are  the  same, 
d 

i.  e.  ,  u  and  thus  x  .  are  the  same,  if 
s  sO 

G(0,t)  -  G(x1,t)  =  3  (y  +  1)  sin  (2t)  (101) 

Referring  to  the  definition  of  G  in  Eq.  (4c),  one  sees  that  there  can  be 
no  plunging  motion  of  the  walls,  but  that  there  is  symmetric  rotation; 
i.e.  ,  G  =  (x  -  xc)p(t)  -  Gu*  From  Eq.  (101),  then,  one  finds  the 

equation  for  p(t),  with  the  result  that 

(x  -  x  ) 

G  =  - - -  3  (Y+l)  sin(2t)  (102) 

X1 

Hence,  if  the  walls  oscillate  with  this  function  G,  if  G^  =  0  so  that  the 
back  pressure  is  constant,  and  if  all  other  conditions  are  as  in  Eq.  (99), 
the  shock-wave  velocity  and  position  will  be  as  shown  in  Fig.  5. 

From  Eqs.  (44)  and  (45),  and  (102),  it  is  seen  that  the  values  of 
c,  (t)  +  G(x,t)  in  h,  are,  upstream  and  downstream  of  the  shock  wave 

Lt  wX 

respectively, 

c,  +  G(x,t)  =  -  —  3  (Y+l)  sin(Zt)  (103a) 

2u  Xj 

(s)  (xl"x) 

cOJ  +  G  (x ,  t )  =  c,Y+  — - -  3  (Y+l)  sin  (2t)  (103b) 

uQ  mu  x  2 

Hence,  the  isotachs  and  pressure  distributions  will  vary  from  those 
shown  in  Figs.  4a,b,c  by  the  time-dependent  terms  in  Eqs .  (103), 
which  occur  in  u  and  p  ,  the  terms  of  order  €^. 
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Numerical  computations  were  carried  out  with  Eqs.  (103)  for 
c^  +  G  in  and  with  all  remaining  conditions  as  in  Eq.  (99).  In  Figs. 
6a,b,c,  isotachs  and  pressure  distributions  downstream  of  the  shock 
wave  are  shown  at  the  same  times  as  in  Figs.  4a,b,c,  for  comparison. 
In  Ref.  [15],  the  flow  upstream  of  the  shock  wave  was  steady  so  no 
isotachs  were  shown  in  that  region;  here  isotachs  upstream  of  the  shock 
wave  are  not  shown  only  because  there  is  nothing  with  which  they  can  be 
compared.  Upon  comparison  of  Figs .  4a,b,c  and  6a,b,c,  one  sees 
that  even  though  the  position  of  the  shock  wave  is  the  same,  the  isotachs 
and  pressure  distributions  are  quite  different.  Two  points  are  of  partic¬ 
ular  importance.  First,  attempts  were  made  to  find  isotachs  with  the 
same  values  as  those  shown  in  Figs.  4a,b,c  so  that  the  differences  in 
shape  could  be  noted,  but  this  was  not  always  possible  because  velocities 
of  the  desired  value  did  not  occur;  hence,  significant  differences  between 
the  two  cases  do  exist.  Second,  it  may  be  noted  in  Figs.  6a,b,c  that 
the  location  of  the  sonic  line  (P  =  1)  does  not  coincide  with  the  throat 
in  steady  flow,  x  =  0,  whereas  it  does  in  Figs.  4a,b,c.  The  reason  for 
this  is,  of  course,  that  the  computations  shown  in  Figs.  6a,b,c  are 
for  the  case  where  the  walls  are  oscillating  so  that  the  position  of  the 
minimum  area,  and  thus  the  sonic  line,  changes  with  time;  this  does 
not  occur  when  unsteadiness  is  caused  by  variations  in  back  pressure. 

It  should  be  noted  that  the  distances  to  the  sonic  line  found  from  com¬ 
putations  of  the  pressure  distribution  and  shown  in  Fig.  6a  were  checked 
using  Eq.  (72),  with  excellent  agreement. 


In  previous  work  described  in  the  references,  the  motion  of  the 
shock  wave  through  the  subsonic  flow  upstream  of  the  throat  was  not 
considered;  hence  numerical  computations  for  u  and  x  stopped  at 
the  throat.  Since  this  analysis  has  been  completed  now,  with  results 
given  in  Eqs.  (97  ),  it  is  possible  to  extend  the  numerical  work  pre¬ 
sented  in  Refs.  [12]  and  [l5]  for  two  example  flows.  Again,  the  results 
are  valid  for  oscillations  either  in  back  pressure,  with  as  in  Eq.  (99) 
or  in  wall  shape  as  in  Eq.  (102).  The  duct  considered  is  that  for  which 
the  wall  shape  is  symmetric  about  the  throat,  so  that  f  (-x)  =  f(x)  with 
f(x)  given  in  Eq.  (99).  The  two  cases  chosen  are  cases  (a)  and  (c)  of 

Fig.  15  in  Ref.  [l 5 3  ,  and  the  results  e  ,r  x  >  0  are  taken  from  this 

su 

report.  For  each  case,  e  =  0.  1,  r  =  150  and  if  oscillating  back  pres- 

(s ) 

sures  are  considered,  G,  =  (Y  +  1 )  4.  5  sin  (2t).  For  case  (a)  x  '  =  0.75 

d  sO 

(g  \ 

and  for  case  (b)  xg^  =  1.5.  The  composite  solution  for  x^  <  0,  Eq. 
(97b),  was  integrated  numerically  using  as  initial  points  those  given  in 
Ref.  [l5]  (Fig.  15)  at  x^^  =  0.  The  results  are  shown  in  Fig.  7.  In 
case  (a),  the  shock  wave  moves  upstream  of  the  throat  very  rapidly, 
moves  out  of  the  duct  and  disappears.  At  a  later  time,  it  forms  again 
at  the  throat,  moves  downstream  and  then  upstream,  again  gaining 
velocity  rapidly  as  it  moves  in  the  subsonic  flow  field.  This  process 
repeats  itself  as  the  back  pressure  (or  walls)  oscillate.  The  shock 
wave  can  never  reverse  its  direction  in  the  subsonic  flow  because  the 
velocity  relative  to  the  wave  must  be  supersonic.  In  case  (c)  a  similar 
situation  occurs  as  the  wave  moves  upstream  and  out  of  the  duct.  As 


shown  in  Ref.  [15],  however,  conditions  are  such  that  after  the  wave 
forms  again  at  the  throat  and  moves  downstream,  it  never  again  moves 
upstream  as  far  as  the  throat;  this  behavior  is  not  shown  here. 

Problem  (2):  or  =  O  (e  "* ) ,  t  =  O  (£  * ) 

In  this  problem,  the  amplitude  of  the  oscillation  of  the  walls 

remains  of  the  same  order  as  in  problem  (1),  but  the  frequency  increases 

such  that  r  becomes  0(£  *);  oscillations  in  back  pressure  again  occur  in 
2 

terms  0(e  ).  For  this  case,  then,  we  set 

a  =  (104a) 

T  =  (kef1  (104b) 

The  equations  for  the  wall  shapes  are  those  given  in  Eqs.  (4). 

Channel-Flow  Solutions;  Velocity  and  Position  of  the  Shock  Wave 

Because  the  wall  shapes  and  thus  the  boundary  conditions  are 
unchanged  from  problem  (1),  the  result  that  the  flow  is  steady  in  first 

order  holds  here  also.  Therefore,  the  shock-wave  velocity  is  again 

2  - 1 
0(e  ),  and  from  Eq.  (7b)  it  is  seen  that,  if  t  =  0(e  ),  x  is  constant 

s  u 

and  the  shock-wave  motion  now  has  an  amplitude  O(e);  thus, 

2 

u  =  kc  dx  ,/dt  +  •  •  •  .  Also,  because  it  is  only  in  second  order  that 
s  si 

a  dependence  upon  time  is  found,  we  note  that,  using  Eq.  (8), 

2 

*T  +2t+‘"  ’°k>> 
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so  that  upstream  of  the  shock  wave  the  Bernoulli  Eq.  (10)  is  again  such 

that  F(T)  =  ^  (Y  +  1 )/ (Y  -  1 )  +  0(€2).  Using  Eq.  (24),  one  can  see  that 

3  1  3 

since  u  (u ,  -  u  )  =  0(e  ),  H  =  H  =  —  (Y  +  1 )/ (Y  -  1 )  +  0(e  );  since  only 
s  d  u  d  u  c 

Z 

terms  to  0(e  )  are  required,  the  stagnation  enthalpy  may  be  considered 
constant  throughout  the  flow  field. 

As  mentioned  previously,  Richey  [  1 8 j  considered  the  problem  of 

flow  in  an  asymmetric  channel  with  oscillations  in  back  pressure  only, 

2  -1  3 

in  the  terms  of  order  e  ,  with  t  =  0(e  ).  Because  or  =  0(e  ),  the 


oscillating  walls  enter  the  boundary  conditions  first  for  and  thus 

affect  only  h^.  The  remaining  terms  in  the  solutions  in  the  channel- 

flow  region  and  the  region  enclosing  the  shock  wave  have  the  same  form 

as  those  found  by  Richey  [is]  .  As  noted  previously,  these  solutions 

have  the  same  form  as  those  found  for  problem  (1),  even  through  r  has 

a  different  order  for  each  case,  because  in  both  problems  the  flow  is 

steady  in  first  order  and  the  walls  are  stationary  to  second  order.  The 

common  solutions  and  therefore  those  that  hold  for  this  problem  are  the 

solutions  for  (t>^  and  <j)^,  contained  in  Eqs.  (37)  and  (39).  Because  only 

2 

these  outer  terms  and  boundary  conditions  up  to  terms  of  order  e  are 
involved  in  the  analysis  of  the  inner  region  behind  the  shock  wave,  the 
solutions  in  that  region  are  also  unchanged  from  those  presented  in 
problem  (1).  Hence,  the  composite  solutions  given  in  Eqs.  (77)  hold 


here  as  well;  moreover,  the  equation  for  u  also  holds,  but  it  represents 


a  differential  equation  for  x^;  i.e.  »  Eq.  (77f)  becomes 


(105) 


4  4k  si  _0_  2 

(Y+l)£2  "  (V  +  1 )  dt  "  (h2x>d  +  (h2x>u-  3  '  Ulu 

2 

It  remains  only  to  find  h^  to  complete  the  solutions  to  0(e  ). 

If  the  expansion,  Eq.  (8),  for  4>^  and  Eqs.  (104)  for  a  and  t  are 

substituted  into  Eq.  (22),  with  H  =  H  =  4-  CY  +  1)/(Y  -  1),  and  into  the 
*  00  £ 

boundary  conditions,  Eq.  (31),  one  finds  the  following  governing  equa¬ 
tion  and  boundary  conditions  for  (j)^,  where  substitutions  have  been  made 
for  the  various  derivatives  of  ^ : 


*3yy  =  2kh2xt  +  <Y+  [*lx  <‘"f"  +  «i  *  +  h2x> 

♦ *21 


♦3y(*.i.«)  =  .*lj(r  +cra 

4,3y(x,-i,t)  =  *lit£;-G(x 


(106a) 

(106b) 

(106c) 


f 


Next,  if  Eq.  (106a)  is  integrated,  the  boundary  conditions  are  applied, 
and  one  of  the  resulting  equations  is  subtracted  from  the  other,  one 
obtains  the  governing  equation  for  h^: 


2khz*.+  «  +  d[+u  £  -  +,x+  -iTTnl"07' 


where  f  and  G  are  defined  in  Eqs.  (37b)  and  (40b),  and  where 
(Y  +  1 )  (|>jx  =  -  f  has  been  used  in  deriving  Eq.  (107).  Upon  com¬ 

paring  this  equation  with  Eq.  (40a),  the  expression  defining  h^  in 
problem  (1),  one  sees  that  the  effect  of  t  being  smaller  by  0(e)  is 


to  introduce  the  time  derivative  of  h^x  into  the  governing  equation. 

It  is  possible  to  separate  the  effects  of  oscillations  in  back 
pressure  and  wall  shape  by  making  the  following  substitution 


‘•’lx  h2x 


lx 


Y+l 


■1 


lx 


where  R(x,t)  contains  the  effects  of  the  wall  oscillations  and  the 
functional  form  of  M  depends  on  the  back  pressure  at  x  =  x^,  as  well 
as  upon  R(Xj,t).  That  is,  if  Eq.  (108)  is  used  for  h^x  and  the  condition 
on  pressure  is  applied  at  x  =  x^,  then,  as  in  problem  (1), 


M(t)  +  R(Xj  ,t)  =  +  Gd(t) 


(109) 


If  Eq.  (108)  is  used  for  h^x  in  Eq.  -(107),  one  finds  the  following 
governing  equation  for  R(x,t) 


fruj  Rt  +  *lx\  =  *lxGx  =  ^lx^ 


Pu  +  P| 

P  =  —  r-—  =  Pft) 


(1 10a) 


(110b) 


where  p^  and  p^,  defined  in  Eq.  (4c),  are  associated  with  the  pitching 
motion  of  the  walls.  Thus,  it  is  seen  that  the  equation  for  R(x,t)  along 
characteristics  and  the  equation  defining  the  one  family  of  character¬ 
istics  are,  respectively, 

A  TJ  (ViH 

(111a) 


*  OLtUa  p 

dt  2k  Vlx  p 

dt  2k  *lx 


(111b) 
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Since  it  does  not  appear  possible  to  find  analytical  solutions  for 
general  wall  shapes,  numerical  computations  must  be  employed.  How¬ 
ever,  for  special  wall  shapes,  it  is  possible  to  find  analytical  solutions 
which  illustrate  the  important  features  of  the  flow  field. 

Analytical  Solutions  for  a  Simple  Case 

Analytical  solutions  for  special  cases,  while  lacking  generality, 
are  valuable  in  giving  a  qualitatively  correct  indication  of  the  physical 
effects  to  be  expected  in  general,  and  in  showing  the  manner  in  which 
these  effects  will  change  as  the  important  parameters  are  varied.  Such 
solutions  can  provide  an  important  guide  for  numerical  calculations, 
serving  as  reference  solutions  exhibiting  the  trends  which  should  be 
anticipated.  In  the  case  considered  here,  the  flow  accelerates  through 
sonic  speed  at  the  throat,  and  a  shock  wave  is  present  further  down¬ 
stream.  The  back  pressure  is  held  constant,  but  the  walls  undergo  a 
simple  harmonic  symmetric  pitching  motion.  The  replacement 
p(t)  =  Pq  e1*  is  made  in  Eq.  (4c),  which  becomes 

G«./-,x-xc)Poei'  <112) 

where  t  =  €  k  T,  so  that  the  (scaled)  frequency  implicit  in  Eq.  (112)  is  k. 
This  form  of  time  dependence  is  particularly  significant  because  an 
arbitrary  periodic  function  of  time  can,  of  course,  be  represented  by 
a  Fourier  series,  and  the  problem  formulated  above  is  linear  with 
respect  to  the  wall  motion,  so  that  solutions  for  more  general  periodic 
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functions  can  be  constructed  by  superposition  of  solutions  having  the 

form  derived  here.  To  simplify  the  problem  as  far  as  possible,  the 

walls  are  taken  to  be  symmetric,  with  the  parabolic  shape  f (0 )  -  f(x) 

1  2 

=  T  (Y  +  1)  (u.  x/x. )  for  -x,  <  x  <  x,  as  given  by  Eq.  (93),  where 
2  loo  1  1  1 

1  2 

f(0)  =  —  (Y  +  1)  u  by  the  choking  condition  (74a);  here  we  have  taken 
Z  I  oo 

Xq  =  -  Xj.  This  shape  might  represent  the  first  term  in  a  series  of 
polynomials  representing  a  more  general  wall  shape.  Analytical  solu¬ 
tions  may,  however,  be  impractical  for  more  complicated  shapes, 
especially  because  the  shape  function  f(x)  enters  the  solutions  in  a  non¬ 
linear  way. 

For  the  case  considered,  the  channel-flow  solution  for  4>^  is 

1  2 

given  by  Eq.  (37a),  with  c.  =  f (0 )  =  —  OY  +  1)  u.  from  Eqs.  (38)  and 

1  c  1  00 

(48)  or  (74a).  The  differential  equation  for  h^  is  then  found  by  sub¬ 
stitution  in  Eq.  (107).  The  results  are,  for  -x^  <  x  <  x^, 

<t>lx=±|uloo,X/Xl  (U3a) 

2kh2xt±(Y+1)(|ulco!/xl)(xh2x)x 

=  P0eiti.(V+l)(!uloo|/x1)3{--£  (Y+l)  +  (|-y)x2} 

(113b) 

where  the  upper  signs  are  used  everywhere  upstream  of  the  shock  wave 
and  the  lower  signs  are  used  downstream.  The  solution  ahead  of  the 
shock  wave,  required  to  be  finite  at  the  throat  x  =  0,  is  easily  found  by 
considering  separately  the  steady-state  part  and  the  simple  harmonic 


unsteady  part.  The  result  is 


l  ui  2  ? 

h2x=6(x)  {“{V  +  1)  +  (3-2Y)xZ} 


P0X1 


(Y+l)  |uloQ|  (l+v2)1/2 


i(t-X) 


(114a) 


v  =  tanX 


2  k  Xj 

(V+l)  |u 


(114b) 


If  a  is  evaluated  from  Eq.  (21b),  it  is  seen  that  the  speed  |u-a|  at  which 

small  disturbances  travel  upstream  becomes  u  -a  =  — (Y+l)eu,  +••• 

oo  oo  L  loo 

outside  the  channel.  The  parameter  v  is  a  reduced  frequency  based  on 

this  velocity  and  on  the  frequency  ke  and  the  length  Xj.  The  phase  lag 

angle  X  increases  from  zero  to  it/ 2  as  v  increases  from  zero  to  infinity, 

1  2  2 

and  X  increases  as  the  magnitude  of  the  wall  curvature  —  (Y+l)u,  /x. 

2  loo  1 

decreases. 

As  for  problem  (1),  a  local  leading-edge  solution  is  required  in  a 

1/2 

thin  vertical  strip  containing  the  edges  and  having  width  0(e  ).  In 

the  present  example,  however,  the  edges  are  wedge-shaped  rather  than 
cusped,  and  so  a  different  solution  is  needed.  The  series  representation 
for  the  perturbation  potential  now  has  the  form 


=  (u  -  l)x  +  €2  4>*  (x^,  y*,  t)  + 

00  J 


x  +  x, 

If  1 


X  = 


(Ke) 


1/2 


y  =  y 


(115a) 


(115b) 
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•jl# 

The  differential  equation  and  boundary  conditions  for  4>^  are 


*>,  +  4>,  =0 

j  *  *  j  *  # 

-XX  ry  y 


(116a) 


4^  (x*,±  l.t)  =Tf,(-x1) 


(116b) 


The  perturbation  velocity  must  decrease  to  zero  as  x*  +  y*  —  oo  out¬ 
side  the  channel  and  must  match  with  the  channel  solution  as  x  r  -*■  oo 
with  -  1  <  y*  <  1.  The  conformal  transformation  (55)  again  maps  the 
flow  onto  the  upper  half  of  the  complex  4  -plane.  The  solution  is 


2  2 

(Ke)1/2(u-u  )-iv  =  ~  f'(-x  )  In  (1  -  -7)  +  ~r—  + 
00  w  1  r  2(4-1) 


(117) 


where  A  is  to  be  determined.  As  z  =  x  +iy  **00  outside  the  channel, 

the  right-hand  side  becomes  e  {A  -  2  f '(-x^  )/ir}/2  u  z  +  •••,  the  form 

for  a  source  at  z  =0,  and  at  the  edges  z  =  ±  i  the  solution  has  an 

inverse  square-root  singularity.  Each  of  these  effects,  if  nonzero,  is 

-1/2 

stronger  by  a  factor  0(6  )  than  for  the  previous  case  of  cusped  edges. 


As  z'"’  *♦  00  inside  the  channel. 


(K€)1/2  (u-u  )  -  i  v  =  €2  f*  ( 
00 


+  -  +  iy 

7T 


.(  X+Xl 

-xiM - 17 

( (K e)  ' 

.  )  12 

■iyj  -  j  € 


e2  A  + 


(118) 


For  the  choked  condition,  the  channel  solution  for  u  =  1  +  4>  will  have 

x 

3/2 

no  term  of  order  e  ,  as  noted  below.  It  follows  that  the  matching 
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with  the  channel  solution  requires  that  A  have  the  constant  value 

2 

A  =  -  f  (-x.)  =  -  -  —  (119) 

it  1  it 

That  is,  A  <  0  and  the  solution  (117)  evaluated  as  z*  *♦  ±i  then  shows 
that  in  this  case  there  is  an  inward  flow  around  the  leading  edges. 

The  source  strength  A  -  2f*  (-x^)/tt  is  zero;  i.  e.  ,  to  the  order  con¬ 
sidered  here,  no  source  term  would  appear  in  the  solution  for  the  outer 
flow.  Thus  the  strength  of  the  source  seen  by  the  outer  flow  has  the 
same  order  of  magnitude  for  both  wedge-shaped  and  cusped  leading 

edges.  This  is  as  expected,  since  the  source  strerigth-sfecuild. npt _  ^ 

depend  on  the  edge  details  but  on  conditions  further  downstream.  Again 
the  channel  solution  has  been  completed  without  knowledge  of  the 
leading-edge  ablution,  but  near  the  channel  entrance  the  flow  is  two- 
dimensional  in  the  first  approximation,  and  the  leading-edge  solution 

is  essential  for  a  correct  flow  description  there.  If  desired,  the  more 

5/2 

complicated  term  of  order  €  in  <f  could  also  be  calculated,  in  the 
manner  of  Ref.  [l £>]. 

In  a  thin  region  near  the  throat  where  x  =  O(e),  it  can  again  be 
shown  that  a  polynomial  solution  in  x  and  y  gives  a  correct  flow  descrip¬ 
tion.  As  in  the  discussion  preceding  Eqs.  (66)  and  (67),  it  is  anticipated 
that  a  nonlinear  differential  equation  may  be  needed,  and  the  wall  bound¬ 
ary  condition  must  also  be  satisfied.  These  considerations  suggest  the 
expansions 
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(120a) 


66 


introduced,  the  differential  equations  combined  with  the  boundary  and 
matching  conditions  would  have  shown  that  the  terms  must  be  identi¬ 
cally  zero. 

The  channel  solution  downstream  of  the  shock  wave  requires 
additional  terms  from  the  solution  to  the  homogeneous  equation  so  that 
the  downstream  boundary  condition  can  be  satisfied.  The  result  is 


2x 


{-  (Y+l)  +  (3-2V)  x*}  - 


(s) 

c2d  X1 


{Y+1)  KJ* 

P0X1  (.  X1  .  x  vi  i (t  +  X ) 

- -  <  1  -  —  exp  iv  In — )>  e  (124) 

u,  |(l+v2)1/2)  X  X1; 


(V+l)  |u.  |(1 

X  CO 

/g  ) 

The  constant  c ^  has  the  same  meaning  as  in  Eq.  (43),  with  a  value 

determined  by  the  constant  prescribed  back  pressure.  If  the  back  pres¬ 
it 

sure  were  to  vary  with  time  as  e  ,  the  integration  constant  chosen  for 
the  solution  to  the  homogeneous  equation  would  be  different.  It  can  be 
seen  that  this  solution  to  the  homogeneous  equation  corresponds  to  the 
function  M  in  Eq.  (108),  having  constant  phase  t  +  X  +  In  (x/Xj)  along 
the  characteristics  dx/dt  =  -  x/v  and  hence  describing  disturbances 
which  travel  upstream  at  this  speed  from  the  channel  exit.  Again  the 
particular  solution  gives  the  direct  response  to  the  motion  of  the  walls. 
The  in-phase  part  has  a  different  sign  than  in  the  upstream  solution 
because  positive  f}(t)  now  implies  a  velocity  decrease,  instead  of  an 


increase  as  was  true  upstream  of  the  shock  wave;  this  is  more  clear 
in  the  limit  as  v  —  0.  The  out-of-phase  part  retains  the  same  sign  as 


r » 


in  the  upstream  solution,  since  it  occurs  as  a  consequence  of  the  4> 


xt 


term  which  has  a  damping  effect  on  the  fluid  motion. 

The  motion  of  the  shock  wave  now  can  be  determined  from  Eq. 
(105).  The  functions  (h?  )  and  (h,  }  are  obtained  by  evaluating  the 

4.X  u  4X 

solutions  (114a)  and  (124),  respectively,  at  the  first  approximation 

x  =  xgg  to  the  shock-wave  position.  The  constant  value  x  is  found  by 

(s ) 

setting  the  time  average  of  dx  ,/dt  equal  to  zero,  so  that  c  '  has  the 

si  id 

same  value  as  in  the  steady-flow  case.  The  results  are 

3 


dx 


si 


‘oxi  e 


(s) 

'2d 

it 


2  V  (V  + 


ii  /|uiJ*soV 


(125a) 


dt 


4|»iJU+v2) 


77 T  i  -  2  i  sin  X  +  — ~  exp  [i  ( X  -  v  In  — — )] } 
1/2  <  Xs0  Xs0  ' 


(125b) 


The  first  term  contains  the  effects  of  the  upstream  and  downstream 

changes  in  flow  velocity  which  result  directly  from  the  wall  motion,  and 

is  180°  out  of  phase  with  the  velocity  ip  (x  -x  )  eXt  of  the  wall  motion. 

0  c 

At  a  value  of  t  when  the  channel  width  is  decreasing  upstream  and  in¬ 
creasing  downstream,  the  flow  velocity  is  increasing  ahead  of  the 
shock  wave  and  decreasing  behind  the  shock  wave,  and  the  first  term 
in  dx^j/dt  implies  that  the  shock  wave  is  moving  in  the  upstream 
direction.  This  term  increases  as  the  reduced  frequency  v  increases, 
since  the  lag  in  the  shock-wave  position  increases  as  v  increases.  The 
second  term  in  dx^/dt  represents  the  effect  of  flow  disturbances  which 


arise  because  the  back  pressure  is  prescribed  and  which  travel  upstream 


at  a  speed  u  -  a  from  the  channel  exit.  The  phase  lag  vln(x,/x  n)  is 

1  sO 

equal  to  the  time  At  required  for  a  disturbance  to  travel  upstream 

from  the  exit  x  =  x,  to  the  shock  wave  x  =  x  _  +  •  •  •  at  a  speed 
1  sO  r 

u  -  a  =  -  (Y+  1)  |uloJ  x/Xj  +  •  •  •  . 

The  simple  nature  of  the  solutions  for  this  special  case  permits  a 

demonstration  of  the  relationships  between  the  solutions  to  problem  (1) 

-2  -1 

and  those  to  problem  (2),  for  which  t  =  0(e  )  and  t  =  0(e  ),  respectively. 

The  general  principle  is  that  the  solutions  to  problem  (1)  evaluated  as 
2 

€  r  -*  0  should  agree  with  the  solutions  to  problem  (2)  evaluated  as  €.t  —  oo  . 

The  present  solutions  can  be  investigated  as  £t  -*  oo  by  taking  the  limit  as 

k  *»  0.  In  this  limit,  the  upstream  solution  (114a)  for  h^  in  fact  has  exactly 

the  form  (40a)  with  c^(t)  +  G(x,t)=  PQx  e^,  which  is  zero  at  the  throat  x  =  0 

so  that  the  condition  (45)  is  satisfied.  Similarly,  as  k  -  0  the  downstream 

solution  (124)  for  h,  has  the  form  (40a)  with  c,  (t)  +  G(x,t)  =  c^J  +  B  (x-x.Je^, 

which  is  constant  at  the  exit  x  =  x^  so  that  the  condition  (44)  is  satisfied. 

Thus  the  solutions  for  the  flow  velocities  in  problem  (1)  are  limiting  cases 

of  the  solutions  for  problem  (2).  This  is  not  true,  however,  for  the 

shock-wave  velocity.  As  k  0,  Eq.  (125b)  shows  that  the  shock-wave 

velocity  u  =  dx  / dT  becomes 
s  s 


u 

s 


(126) 


The  shock-wave  velocity  in  problem  (1),  as  expressed  by  Eqs.  (79),  can 
2 

be  expanded  as  e  t  •*  0  by  noting  that  u^  approaches  a  constant  value 
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because  the  amplitude  of  the  shock-wave  oscillation  becomes  small. 

(s) 

The  value  of  C2^  is  again  given  by  Eq.  (125a),  and  substation  for  G 

shown  that  Eq.  (79a)  again  gives  the  value  (126)  for  u  ,  as  it  must. 

s 

Problem  (3):  a  =  0(e^)t  r  =  0(6  *) 

In  this  problem,  the  amplitude  of  the  oscillations  of  the  wall  is 
of  the  same  order  as  the  thickness  of  the  wall  (blade)  itself,  and  is 
therefore  an  order  of  magnitude  larger  than  in  problems  (1)  and  (2). 
The  frequency  is  the  same  as  in  problem  (2),  i.e.  ,  t  =  0(£_1).  Thus, 
for  this  case,  we  set 

a  =  €2  (127a) 

t  =  (kef1  (127b) 

The  equations  for  the  wall  shapes  are  again  those  given  in  Eqs.  (4). 


Governing  Equations,  Channel- Flow  Region 

Just  as  in  problems  (1)  and  (2),  application  of  the  boundary  con¬ 
ditions  to  the  second-order  terms  leads  to  the  derivation  of  the  govern¬ 
ing  equation  for  the  first-order  term  in  the  velocity  potential  4>j. 

Because  the  boundary  conditions  do  not  involve  first-order  terms, 

2 

4»j  is  again  independent  of  y;  however,  because  terms  of  order  e  in 
y  are  time -dependent,  the  first-order  terms  are  no  longer  steady, 


w 


but  depend  upon  time.  As  a  result,  u  =  O(e),  and  so  from  Eqs.  (76) 

s 


and  (127b)  u  =  k  6  dx  n/dt  + 
s  su 


That  is,  x  _  =  x  „(t)  and  the 
sO  sO 
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magnitude  of  the  oscillations  of  the  shock-wave  position  is  0(1);  large 

excursions  of  the  shock  wave  can  occur.  In  addition,  the  oscillations 

in  back  pressure  can  occur  in  the  term  of  order  € ,  rather  than  the 

2 

term  of  order  €  as  was  the  case  in  problems  (1)  and  (2). 

In  order  to  derive  the  governing  equations  for  the  4>.  in  Eq.  (8), 

which  again  is  used  to  represent  the  perturbation  potential,  it  is 
2 

necessary  to  find  a  ;  that  is,  it  is  necessary  to  evaluate  F(T)  in  Eq. 

(10).  For  this  case  it  is  again  convenient  to  use  the  energy  equation 

written  in  terms  of  the  stagnation  enthalpy,  Eq.  (15),  and  to  expand 

H  as  in  Eq.  (17a)  with  the  governing  equation  for  given  by  Eq.  (17b). 

Comparing  Eqs.  (16)  and  127b),  we  see  that  n  =  1  in  Eqs.  (16)  and 

(17a)  for  this  case.  In  addition,  Eqs.  (13),  (19),  (20),  and  (21)  hold 

2 

here  for  Pj,  P^,  Hj,  F(T),  and  a  ,  respectively.  In  Eqs.  (20) 

2  2 

and  (21a),  and  thus  in  a  since  H  -  H  =  e  H,  ,  the  function  e.  occurs 

oo  1  °1 

and  must  be  evaluated  both  upstream  and  downstream  of  the  shock 
wave.  Although  g^  =  gj(y,t)  as  a  function  of  integration  of  Eq.  (17b), 
it  is  seen  from  Eq.  (21a)  that  it  is  at  most  a  function  of  t.  The  evalua¬ 
tion  of  gj  is  more  easily  accomplished  by  considering  the  integration 
of  Eq.  (17b)  at  a  given  instant,  i.  e.  ,  by  evaluating  a  definite  integral. 
Thus,  upstream  of  the  shock  wave, 

Hj  (x,t)  -  Hj  (-oo,  t)  =  -  k  («j>lt(x,t)  -  4>lt  (-oo,  t) )  . 

But  Hj  (-oo, t)  =  0  by  definition,  and  we  arbitrarily  choose  the  unknown 
function  of  time  in  the  potential  to  be  zero,  so  (-oo,t)  =  0.  Hence, 


1  1  5W 


upstream  of  the  shock. 


Hj  =  -  k  <>lt  (128a) 

gj  =  0  (128b) 

In  order  to  perform  the  same  kind  of  calculation  downstream  of 
the  shock  wave,  it  is  necessary  to  calculate  the  jump  in  Hj  across  the 
wave.  From  Eqs.  (24)  and  (30),  it  is  seen  that  to  the  order  desired, 

H  =  H  +  u  (u,  -  u  )  +  •••  (129) 

d  u  s  d  u' 

where,  again,  the  subscripts  u  and  d  refer  to  positions  immediately 
upstream  and  downstream  of  the  shock  wave,  respectively.  If  Eqs. 

(5a)  and  (17a)  are  used  for  the  expansions  of  u  and  H,  then  £q.  (129) 
may  be  used  to  show  that 

Hld  =  Hlu  +  ^  <uld  '  ul»>  ,130,) 

Hlu--k,Vu  (I30b> 

where,  since  the  flow  is  unsteady  in  first  order,  u  must  be  O(e),  as 

s 

mentioned  previously.  If  now  Eq.  (17b)  is  integrated  as  a  definite 
integral,  at  a  given  instant,  from  the  shock  wave  to  a  point  somewhere 
downstream,  it  is  found  that  downstream  of  the  shock  wave 


Hi 3  Hid  - k  <ht  -  <Vd» 


*1  1  Hld  +  k  (Vd 


where  will  be  calculated  later. 


(131b) 


If  Eq.  (21b),  with  H  -  H  =  e^H  ,  and  Eq.  (8)  are  substituted 
into  Eq.  (9),  with  T  =  rt  and  t  as  in  Eq.  (127b),  one  obtains  the 
following  governing  equations  for  4>^,  4>2,  ■ 


(132a) 

(132b) 


hyy 

*2yy  *  2k*l*»  +  (Y+  »  hx  +  lx*  +  2  hy+lxy 
*3yy  =  ^  *ltt  +  2  »  *2xt  -  (V  -  1)  H,  +  2k 


+  +  <V-1»*lx4.2yy+  ^  <£♦,„  (132=) 


The  corresponding  boundary  conditions  are  found  from  Eq.  (4),  with 


a  =  €  ,  and  Eqs.  (5)  and  (31);  thus 


+2y(r»1*t,=  ’(fi‘G«) 


4»2y(x, -i,t)  =  f;  G|x 

/v  1  f  \  -  -  ^ 

lx  U  UX 
(j>3y  (x,  t)  =  «|>lx  (fff  -  G|x)  -  k  Gft 


4>3y  (X»  !*  t]  =  •  +lx  (fu  "  Gux)  +  kGut 


(133a) 

(133b) 

(133c) 

(133d) 


where  <t>.  -  u.  and  4>.  -  v. . 

Tix  x  ly  l 


Since  4>,  =  0  at  either  boundary,  then,  from  Eq.  (132a), 

ly 


a  =  <j>,  (x,t)  =  0  and  so 

ly  iy 


=  <t>j  C*»t) 


With  this  result  for  <j>  ,  Eq.  (132b)  may  be  integrated  to  give 

2 

4,  =  ((Y+i)  4>lx <t>lxx  +  2 k 4>lxt)  *2  +  s2y  +  h2(x,t) 


(134) 


(135) 
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After  applying  the  boundary  conditions  on  4>  given  by  Eqs.  (133a)  and 
(133c),  one  obtains  two  equations  which  may  be  combined  to  give 


g2(x..)=i[£-  (136a> 

Zk*b»t(,tl)^*ta=-i’4Gx  <136bl 


where  f  and  G  are  defined  in  Eqs.  (37b)  and  (40b),  respectively.  Eq. 
(136b)  is  seen  to  be  the  governing  equation  for  <j>^ ;  as  mentioned 
earlier,  this  case  is  different  from  those  considered  in  problems  (1) 
and  (2)  because  here  <j>^  depends  upon  time.  It  may  be  noted  that 
whereas  the  effects  of  the  oscillating  walls  are  felt  in  the  governing 
differential  equation  itself  through  the  function  G,  the  oscillating  back 
pressure  is  introduced  either  through  initial  or  boundary  conditions. 

In  order  to  complete  the  solutions  to  second  order,  it  is  neces¬ 
sary  to  find  the  governing  equation  for  h^x.t).  As  in  the  previous 
cases  considered,  this  equation  is  found  by  integrating  Eq.  (132c)  once, 
applying  boundary  conditions  given  by  Eqs.  (133b)  and  (1 33d),  and  sub¬ 
tracting  one  of  the  resulting  equations  from  the  other.  The  resulting 
equation  is,  after  some  rearranging, 


2kh2xt  4  (V+  *)  (*lxh2x)x  1  "  htt  '  t  Gxxt  4 

-  2k*i*  hxt  hxx  -  *i*,f-G)*+kGt  <137) 
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Equation  (137)  may  be  simplified  somewhat  by  making  substitutions 


for  and  <t>2xt*  Thus,  if  Eq.  (136b)  is  integrated  over  x,  again  in 
terms  of  a  definite  integral  with  (-00, t)  =  0  as  mentioned  above, 
and  the  resulting  equation  is  differentiated  with  respect  to  t,  one  finds 


that 


k  <|> 


ltt 


(V  +  1) 
2 


^lx  ^Ixt 


(138) 


where  (-00, t)  =  u,  =  constant  here  and  where  F  =  G  =  0  upstream 

lx  loo 

of  the  channel,  where  x  <  x^.  If  Eq.  (138)  is  substituted  into  Eq.  (137), 
and  if  Eq.  (136b)  is  used  for  $  in  the  resulting  equation,  the  gov¬ 
erning  equation  for  h_  is  found  to  be,  finally. 


fH*,2*,  ♦.  «-G> 

4  lx  lxx  4  lx  x 


CL- 


(139) 


The  last  equation  needed  for  solution  of  this  problem  to  the  order 
desired  is  that  for  the  velocity  and  position  of  the  shock  wave.  As  in 
the  problems  considered  previously,  this  relation  is  found  from  the 
jump  conditions  which  hold  across  the  moving  shock  wave.  The  rele¬ 
vant  equations  are  (30d)  and  (30e),  where  u  and  H  on  the  right-hand 
side  of  Eq.  (30e)  may  be  evaluated  either  both  upstream  or  both  down¬ 
stream  of  the  shock  wave,  as  a  result  of  Eq.  (24);  here  conditions  up¬ 
stream  of  the  wave  are  used.  With  Eqs.  (5a)  for  u,  (7b)  for  ug ,  and 
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A 


/dt  and 


(140a) 


(140b) 


Thus,  it  is  seen  that  to  first  order  the  jump  conditions  across  the 

shock  wave  may  be  satisfied  as  long  as  Eq.  (140a)  is  satisfied.  The 

only  point  at  issue  is  the  initial  condition  in  that  x^q  must  be  specified 

at  a  given  time;  it  is  not  possible,  for  example,  to  use  Eq.  (140a)  to 

obtain  the  steady- state  value  of  x  n,  when  dx  n/dt  =  0.  On  the  other 

sO  sO 

hand,  it  is  clear  that  Eq.  (140b)  cannot  be  satisfied  in  that  x  =  x  .  (t) 

si  si 

whereas  u^  and  u^u  both  depend  upon  y  as  well  as  upon  time.  Thus, 

the  jump  conditions  cannot  be  satisfied  to  second  order  by  the  channel- 

flow  solutions.  As  in  previous  cases,  it  will  be  necessary  to  consider 

an  inner  region  about  the  shock  wave  which  is  now  moving  at  u  =  O(e). 

s 

In  this  region,  solutions  will  have  to  satisfy  the  local  jump  conditions 
across  the  shock  wave  and  match  with  the  outer  channel-flow  solutions. 
The  inner  region  has  not  been  analyzed  as  yet  for  this  case.  Hence  the 
problem  formulation  described  here  allows  solution  only  to  first  order, 
using  Eqs.  (136b)  and  (140a). 


76 


Simplified  Formulation  for  a  Special  Case 

Simplifications  leading  to  an  ordinary  differential  equation  for 
the  shock-wave  motion  can  be  achieved  in  a  special  case  analogous  to 
that  considered  for  Problem  2.  Although  for  the  present  problem  a 
complete  analytical  solution  can  not  be  obtained,  the  example  appears 
to  provide  a  useful  test  case  for  numerical  calculations,  to  show  the 
general  effects  of  changes  in  the  parameters  with  numerical  integra¬ 
tion  needed  only  for  the  shock- wave  motion  and  not  for  the  entire  flow 
field. 

In  this  special  case  the  flow  accelerates  through  sonic  speed,  at 
a  location  which  in  general  no  longer  coincides  with  the  location  of  the 
instantaneous  minimum  cross-section  area,  and  a  shock  wave  is  pres¬ 
ent  further  downstream.  The  walls  have  the  symmetric  parabolic 
shape  given  by  Eq.  (93)  and  are  given  a  simple  harmonic  pitching  os¬ 
cillation  as  shown  by  Eq.  (112),  where  t  =  e  k  T  as  before,  so  that  the 

(scaled)  frequency  is  again  k.  The  difference  is  that  the  amplitude  of 

2  3 

the  wall  oscillations  is  0(£  )  rather  than  0(6  );  it  follows  that  the  sonic 
line  and  the  shock  wave  undergo  displacements  which  are  0(1)  rather 
than  O(e). 

With  the  choices  (93)  and  (112)  for  the  functions  f(x)  and  G(x,t), 
the  differential  equation  (136b)  can  be  rewritten,  with  4>jx  =  Uj  for 
simpler  notation,  as 

2ku  +  (Y+  1)  u,  u  =  (Y+  l)(u  /x  )2x  +  B  eAt  (141) 

It  1  lx  loo  1  0 
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For  an  arbitrary  right-hand  side,  this  equation  can  be  integrated 
numerically  along  characteristics  dx/dt  =  ■-  (Y  +  1 )  u  /k.  In  the  pres- 

Lt  x, 

ent  case,  in  order  to  obtain  an  analytical  result,  a  coordinate  trans¬ 
formation  from  x,t  to  r,s  is  carried  out,  such  that  s  is  measured  along 
characteristics  and  so  r  is  constant  along  characteristics;  i.  e.  ,  the 
value  of  r  identifies  a  characteristic.  If  we  take  dt/ds  =  1,  the  differ¬ 
ential  equation  ( 141)  is  replaced  by  the  system 


at, 

8s 


1 


(142a) 


8x  _  Y+  1 
8s  =  2k  U1 


8s 


y  +  i 

2k 


u,  “ 
*1 


it 

e 


(142b) 


(143c) 


It  is  convenient  to  take  s  =  0  at  t  =  0,  so  that  Eq.  (142a)  gives  simply 
t  =  s.  The  general  solutions  for  x  and  u^  can  be  written  as 


x  =  A(r)  e 


(s-r  )/v 


+  B(r)  e 


-  (s-r)/v 


P0*l 


IS 


(Y+Du^d+v2) 


(144a) 


u 


1 


{A(r) 


(s-r)/v 

e 


B  (r ) 


e-(s-r)/v} 


1  P0X1  V 


is 


(Y+  1)  lujJ  (1  +v  ) 


(144b) 


where  A(r)  and  B(r)  are  to  be  determined  from  initial  conditions  for 
x(r,s)  and  u ^ (r ,  s ),  and  v  is  defined  by  Eq,  (314b).  For  the  flow 
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downstream  of  the  shock  wave,  dx/dt  <  0  along  characteristics,  so  that 


disturbances  propagate  upstream  from  the  exit  cross-section  x  =  x 


1 


toward  the  shock  wave,  and  initial  conditions  should  then  be  specified 
at  x  =  Xj,  We  can  choose  r  =  s  at  x  =  x^ ,  so  that  the  initial  conditions 


are 


Uj  (r ,  r )  =  uQ  , 


x  (r,r)  =  x. 


(145a, b) 


where  u^  is  determined  Ky  the  downstream  pressure;  we  assume  that 
the  0(c)  term  in  the  pressure  at  x  =  x^  can  be  specified  arbitrarily. 
The  solutions  for  x  (r,s)  and  u^  (r,s)  downstream  of  the  shock  wave 
are  then 


-  i+r 

X,  2  1  u 


Up  V0<l+i>>  ir. 

I  +  2  2  6  I  e 

loo1  CY+  1)  u.  (1  +v  )  ; 

1  00 


2«  In. 


X,  (1  -  IV)  .  .  .  .  . 

1  0  lrl  -(s-r)/v 

+  - X - T-  e  >  e 


loo'  (Y+l)u,2  (l+vZ) 
loo 


*1P0 


is 


(Y+l)  (1+v2) 


(146a) 


u 


“1  ll„  °0  .  ir  1  (s-r)/v 

2  KtJ  (V  + 1 1  u 2  (1+v2)  ' 

1  00 


loo 


u„ 


x!  P0  f1  -  iv) 


ir  1  -(s-r)/v 


Koo'  (Y  +  1 )  u  2  (1+v2) 
loo 


e  }e 


lv 


*1  ^C) 


(Y+l)  u  2  (1  +  v2) 
1  00 


e1S  (146b) 
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For  the  flow  ahead  of  the  shock  wave,  the  characteristics 


dx/  dt  =  —  (\  +  1)  Uj  have  positive  slope  at  points  in  the  x,t  plane  where 
Uj  >  0  and  negative  slope  at  points  where  <  0,  so  that  small  dis¬ 
turbances  are  always  propagated  in  a  direction  away  from  the  sonic 
line  where  =  0.  The  first  approximation  to  the  location  of  the  sonic 
line  will  be  denoted  by  x  =  x^(t),  where  x*  is  the  value  of  x  for  which 
Uj  =  0.  We  will  assume  tentatively,  subject  to  a  later  check  for  self- 
consistency,  that  Uj  =  0(x-x^)  as  x  -*■  x#,  just  as  for  steady  flow;  that 
is,  3u,/3x  is  neither  zero  nor  infinite  at  x  =  x*.  It  then  follows  that 
dx/dt  =  0(x  -  x*)  along  characteristics  as  x  -*  x^  and  therefore  also 
s  -*  -oo  as  x  •*  along  a  characteristic.  Thus  necessarily  B(r)  =  0, 
and  the  solution  for  u^  found  from  Eqs.  (144a, b)  is  in  fact  simply  pro¬ 
portional  to  x  -  x^(t): 


_^1 _ x-  *„(*> 


loo 


(147a) 


x,.(t) 

T* 


Pp*!2  “  -lv)  it 

2  Z  6 

(Y+l)  (1  +  vC) 


(147b) 


The  result  of  course  is  consistent  with  the  assumption  made  above. 

The  instantaneous  cross-section  area  of  the  channel  has  a  minimum  at 

the  location  x  =  x  .  where  the  right-hand  side  of  Eq.  (141)  is  zero, 
min 

By  comparison  with  Eq.  (147b)  it  is  seen  that  x*  =  x  .  /(I  +  iv),  and 

min 

so  the  sonic  line  no  longer  occurs  at  the  cross-section  having  minimum 
area.  Since  0  <  arg  (1  +  iv)  <  -tt  /  Z  for  0  <v  <  oo  ,  the  displacement  of 
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the  sonic  line  lags  the  wall  motion  by  a  phase  angle  between  zero  and 
w/2;  the  amplitude  in  Eq.  (147b)  decreases  as  the  frequency  increases, 
with  *  0  as  v  —  oo  . 

A  further  check  concerning  the  flow  near  the  sonic  line  x  =  x^ 

can  be  made  by  observing  that  the  term  2  k  <|>  in  the  differential  equa¬ 
l/2 

tion  becomes  of  the  same  order  as  d>  when  x  -  x.,.  =  0(e  ).  A  local 

yy 

solution  can  be  obtained  in  the  form 


<i>  =  €2  <|>  *  (x*,  y,  t)  + 


*  X  ~  X* 


.1/2 


(148a) 

(149b) 


One  finds  then 


dx 

4>*  +  2k 

2yy  dt 


*  * 

*2x*x*  =  ° 
u2 


4>*  (x*,  +1,  t)  =  ±  (  (Y+  1)  — ~  x^  +  (3q  e11) 
7  X1 


(149a) 


(149b) 


The  solution  which  matches  correctly  as  x  •*  i  oo  with  the  channel 
solution  evaluated  as  x  -*■  x^t)  is 


u. 


«t>2*  (**#  y»  t)  =  j  y2  ( (Y+  1)  x;;.  +  (3q  e11) 

x! 


,  I  lllsl 

2  xj 


x*2  +  h2  (x#,  t)  (150) 


where  the  function  h^x.t)  is  not  yet  known,  but  would  be  found  by  solu¬ 
tion  of  Eq.  (137).  The  result  (150)  simply  says  that  the  largest  terms 
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in  the  channel  solution  remain  correct  near  the  sonic  line  x  =  x^(t)  +  •  •  •  , 
and  that  no  special  solution  is  needed  there. 


The  shock-wave  velocity  dx^/dt  is  given  by  Eq.  (140a)  in  terms 
of  the  flow  velocities  u^  and  u^  immediately  upstream  and  downstream. 

Since  x  is  continuous  at  the  shock  wave,  one  more  condition  is  available. 

To  use  this  condition  we  can  also  define  r  to  be  continuous  across  the 
shock  wave,  so  that  a  characteristic  arriving  at  a  given  instant  from  the 
upstream  side  has  the  same  value  of  r  as  the  characteristic  which 
reaches  the  shock  wave  from  the  downstream  side  at  the  same  instant. 

From  the  continuity  of  x  (r,s),  the  function  A(r)  in  the  upstream  solu¬ 
tions  for  x  and  u^  is  then  determined,  and  it  is  easily  found  that 

! 

U,  -uu  /  un  Pn  xi  ^  -iv)  o\/„ 

-!s - r*  =  -  ‘i — +  - °-  1  ■; - r  elRWR-s)/v  (151) 

Kool  1  lUlooi  (v  +  1 )  U ,  (1+0  ' 

loo 

The  shock-wave  velovity  dxg^/dt  can  be  expressed  in  terms  of  the 
shock -wave  position  R(s)  if  the  solution  for  x(r,s)  is  evaluated  at 
r  =  R(s )  and  then  differentiated: 


9x  9x 

r=R  +  8r 


dR 
r  =  R  ds 


(152) 


With  the  substitution  (142b)  for  x(r,s),  it  is  found  from  the  first  approxi¬ 
mation  (140a)  to  the  shock-polar  equation  that 


Substitution  in  Eq.  (153)  gives,  finally,  an  ordinary  differential  equa¬ 
tion  for  R(s) 


Eq*1  »-iWeiR  1  A  \ 

(Y+l)u2  (1  +  v2))  \dR  / 

loo 


2(s-R)/v 
e  + 


P 

(Y 


0xlelR  fiv(i-i 

+  l)u2  \  1  +  V2 

loo 


vl  Z(s-R) /vi 


(154) 


where  it  is  understood  that  the  real  part  is  to  be  taken.  For  the  sim¬ 
plest  example  the  exit  velocity  is  the  same  as  the  entrance  velocity,  so 

that  u.  =  u,  <  0;  the  remaining  parameters  are  v  and  x,  /[(V+l )  u  2  ]. 
0  loo  0  1  loo 

If  an  initial  condition  is  specified,  numerical  integration  of  Eq.  (154) 

would  give  the  shock-wave  location  in  the  form  r  =  R(s);  substitution  in 

the  solution  (146a)  for  x  (r,s)  then  would  allow  calculation  of  x  (t). 

so 

Presumably  some  sort  of  criterion  for  the  existence  of  periodic  solu¬ 
tions  could  thereby  also  be  established.  As  a  preliminary  step,  quali¬ 
tative  imformation  about  the  solutions  should  be  sought  by  consideration 
of  Emiting  cases  and  by  study  of  the  singular  points  of  the  differential 
equation. 
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HI.  SUMMARY 


The  extension  of  previous  work  on  unsteady  transonic  flow  in  two- 
dimensional  channels  to  include  effects  of  oscillating  walls  has  proven 
to  be  relatively  straightforward,  for  the  cases  considered,  insofar  as 
formulation  of  the  problem  is  concerned.  Thus,  it  has  been  possible  to 
consider  problems  where  the  unsteadiness  may  be  caused  by  any  mix¬ 
ture  of  oscillations  in  wall  position  and  back  pressure.  The  essential 
difference  between  the  two  lies  in  the  fact  that  when  the  walls  are 
stationary,  oscillations  in  back  pressure  cannot  cause  unsteadiness 
upstream  of  the  shock  wave,  whereas  when  oscillations  in  wall  position 
occur,  the  entire  flow  is  time-dependent.  Although  it  may  be  possible 
to  find  equivalent  impressed  oscillations  for  the  walls  and  back  pres¬ 
sure,  such  that  the  same  shock-wave  velocity  and  position  result  (e.g.  , 
problem  1),  the  distributions  in  flow  properties  are  then  not  the  same. 
Hence,  as  a  general  result,  it  is  not  possible  to  find  equivalent  im¬ 
pressed  oscillations  which  give  the  same  instantaneous  flow  fields. 

The  problems  considered  involve  both  small  and  large -amplitude 
shock-wave  oscillations;  in  the  latter  case  motion  of  the  shock  wave 
throughout  the  whole  duct  is  allowed.  In  general,  derivation  of  the 
governing  equations  for  the  flow  properties  has  been  carried  out  to  the 
point  that  solutions  valid  to  second  order  may  be  obtained.  Where 
possible,  analytical  solutions  have  been  presented.  Where  this  was 
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impossible,  analytical  solutions  for  special  cases  have  been  given  and 
the  problems  formulated  for  numerical  calculation. 

Analysis  of  the  thin  region  about  the  leading  edge  has  extended 
previous  analyses  in  that  inlet  problems,  as  well  as  straightforward 
channel  flow,  may  now  be  considered.  Although  exterior  cowl  effects 
were  not  taken  into  account,  because  specific  shapes  would  be  needed, 
their  addition  would  be  straightforward. 

One  of  the  more  interesting  results  found  as  a  result  of  the 
detailed  analysis  of  the  throat  region  is  the  possibility  of  large  excur¬ 
sions  of  the  sonic  line,  and  the  fact  that  the  sonic  line  need  not,  then, 
occur  at  the  instantaneous  position  of  minimum  cross-sectional  area. 

The  application  to  a  duct  with  boundary  layer  which  separates 
downstream  of  a  shock  wave  such  that  the  size  of  the  bubble  (and  thus 
the  effective  channel  wall  shape)  oscillates  with  time  was  considered 
in  detail  only  for  problem  1.  However,  this  illustrated  the  method  of 
solution  well  enough  that  further  demonstrations  were  considered  un¬ 
necessary. 

Detailed  numerical  computations  were  given  only  for  problem  1. 
Only  example  solutions  with  simple  wall  shapes  and  harmonic  oscilla¬ 
tions  were  shown  for  problem  2,  because  more  complicated  shapes  and 
conditions  would  have  necessitated  extensive  numerical  work  which 
would  have  added  little  understanding.  Such  is  not  the  case  for  prob¬ 
lem  3,  where  the  simple  solution  presented  allowed  only  calculation  of 
the  velocity  field  and  not  the  velocity  and  location  of  the  shock  wave; 
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£ 


Moreover,  this  problem  is  important  enough,  involving  unsteadiness 
in  lowest  order  perturbation  as  it  does,  that  it  is  recommended  that 
detailed  and  fairly  extensive  numerical  calculations  should  be  carried 


out  in  future  work. 
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Figure  2.  Coordinate  system  and  notation  for  moving 
shock  wave. 
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Figure  4.  (c)  t  =  2.  35618 


Isotachs  and  wall  and  centerline  pressure  distributions  for  +  G  as  in  Eqs. 
(103)  and  all  remaining  conditions  as  in  Eq.  (99). 

(a)  t  =  0.7854 
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